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Abstract: A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread 
across the celestial sphere are combined, allows investigation of "global" phenomena such as a back- 
ground of gravitational waves or instabilities in atomic timescales that produce correlated timing resid- 
uals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the 
PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond 
pulsars is being observed at three radio-frequency bands, 50cm (~ 700 MHz), 20cm (~ 1400 MHz) and 
10cm (~ 3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced 
in early 2005. This paper describes the systems used for the PPTA observations and data processing, 
including calibration and timing analysis. The strategy behind the choice of pulsars, observing pa- 
rameters and analysis methods is discussed. Results are presented for PPTA data in the three bands 
taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less 
than lfis for the best band after fitting for pulse frequency and its first time derivative. Significant 
"red" timing noise is detected in about half of the sample. We discuss the implications of these results 
on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the 
Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data 
with earlier Parkes timing data for these pulsars. 

Keywords: pulsars: general — instrumentatiommiscellaneous — methods:observational — gravita- 
tional waves 
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1 Introduction 

Pulsars have many intriguing properties, but their most 
important attribute by far is the remarkable stability 
of the basic pulse periodicity. Having these "celes- 
tial clocks" distributed throughout the Galaxy (with 
a few in our nearest neighbour galaxies, the Magel- 
lanic Clouds), many of them members of binary sys- 
tems, makes possible a range of interesting and im- 
portant applications. The best known of these is the 
detection of orbital decay in the original binary pul- 
sar, PSR B1913+16 ijHulse fc Tavlorlll975h . which pro- 
vided the first observational evidence for the existence 
of gravitational waves (GWs) and showed that the 
rate of energy loss was in accordance with the pre- 
dictions of Einstein's gene ral theory of relativity (GR) 
(Wcisbcrg fc Tavlorl 12003 ) . But there are many oth- 
ers. For example, the pulse dispersion due to free elec- 
trons along the path to the pulsar can easily be mea- 
sured and used to study the distribution of ionised 
gas in our Galaxy and, potentially, the intergalactic 
medium. Precise positions, proper motions and even 
the annual parallax of pulsars can be measured using 
pulsar timing. Careful study of the timing of binary 
pulsars has revealed a range of orbital perturbations 
which not only give important information about the 
formation and evolution of the binary systems, but 
also allow sensitive tests of gravitational theories. In 
particular, the discovery and subsequent timing obser- 
vations of the firs t-known double-pulsar system, PSR 
J0737-3039A/B (|Burgav et al.l2003l : lLvne et al.il2004l ). 
has allowed four independent tests of GR and given the 
most pr ecise verification so f ar of GR in the strong-field 
regime (jKramer et al.ll200a ). 

One of the major goals of current astrophysics is 
the direct detection of GWs. Detection and study of 
these waves would open up a new window on the early 
Universe and the physics of extreme gravitational in- 
teractions. Enormous effort is going into the construc- 
tion of systems such as the Laser Int erferometer Gravi- 
tatio nal Wave O bservatory (LIGO) |Abramovici et al.l 
1992) and Virgo l|Acernese et al. 2004) which are sensi- 
tive to GWs with frequencies in the range 10 - 500 Hz. 
Initial versions of these systems are now operating and 
have placed limits on the amplitude of GWs from sev- 
eral types of astrophysical source (e.g., lAbbott et al.l 
2006). Systems with improved sensitivity are being 
developed. Advanced LIGO is due for completion in 
2014. KAGRA, a 3-km underground detector with 
cryogenic ally-cooled optical systems is being constructed 
in Japan (|Somiva fc for the KAGRA Collaborationll201lf ) 
and the Einstein Telescope, a thi rd-generation detec- 
tion system, is under development (|Punturo et al.ll2010f ). 
The Laser Interferometer Space Observatory (LISA) 
has now evolved into eLISAy which will be sensitive 
to GWs with frequencies in the range 0.1 - 100 mHz. 
Finally, space observatories exploring the cosmic mi- 
crowave background such as the Wilkinson Micro wave 
Anisotropy Probe (W MAP) dBennett et al]|2003l ') and 
the Planck Surveyor (jTauber et al.ll2010l ) have as one 
of their goals the detection of the B-mode polarisation 
signature of primordial GWs; these have a local fre- 

1 Scc http:// www.elisa-ngo.org 



quency of ~ 10" 9 Hz (|Crittenden et al.lll995h. So far 
only upper limits have been obtained (e.g. jLarson et ail 

Hqii|). 

Because of the great intrinsic stability of their pul- 
sational periods (P) or pulse frequencies (y = 1/P), 
precision timing observations of pulsars and, in partic- 
ular, millisecond pulsars (MSPs), can in princi ple be 
used to detect GWs p ropagating in our Galaxy (|Sazhinl 
Il978l : lDetweilerl|l979T l . Observed pulse frequencies are 
modulated by GWs passing over the pulsar and the 
Earth; the net effect is the difference between these two 
modulations. Pulsar timing analyses give the differ- 
ences, commonly known as timing residuals, between 
observed pulse times of arrival (ToAs), normally re- 
ferred to the barycentre of the solar system, and the 
pred ictions of a model fo r the pulsar properties (see, 
e.g.. lEdwards e t al. 2006). Pulsar timing is therefore 
essentially a phase measurement and is sensitive to the 
integrated effect of pulse- frequency modulations. Con- 
sequently, it is most sensitive to long-period modula- 
tions out to roughly the data span. This is typically 
many years, corresponding to frequencies in the range 
1-30 nHz. Pulsar timing experiments are therefore 
complementary to the ground-based and space-based 
laser interferometer systems. Since the intrinsic value 
of v and its rate of change, v, are a priori unknown, 
they must be solved for in the timing solution. There- 
fore, any external modulation which affects only these 
parameters is undetectable^ 

Different astrophysical sources are likely to dom- 
inate the GW spectrum in the different bands. For 
example, the GW sources most likely to be detected 
by ground-based interferometer systems are the final 
stages of coalescence of double-neutron-star binary sys- 
tems. For eLISA the most promising source is a simi- 
lar coalescence of super-massive binary black holes in 
the cores of distant galaxies. At an earlier stage of 
their evolution, these same super-massive black-hole 
binary systems generate a stochastic background of 
low-frequency GWs. The expected spectrum of this 
GW b ackground can be described by a power-law re- 
lation (|Jenet et al.ll2006T ): 

hc{f)=As {h^) (1) 

where h c is the characteristic GW strain over a loga- 
rithmic frequency interval at frequency / and /i yr is 
1/lyr. Predictions for the rate of super-massive black- 
hole binary coalescence events following galaxy merg- 
ers result in values of the GW background amplitude 
A s that are pot entially detectable by pulsar timing ex- 
periments (e.g., [Jaife_fc_13ackejj l2003l: IW vithe fc Loebl 
120031 ; lEnoki fc Nagashimal 120071 ; ISesana et alj l200tf ). 
For such a GW b ackground, the exponent a is ex- 
pected to be -2/3 l|Pninnevll200lh . Other potentially 
detectabl e sources include cosmic strings in the early 
universe IjDamour fc Vile nkin 2005: ISiemens et all2007r 
and fluctuations in the primordial Universe (jGrishchukl 

2 An exception to this is where the effect is sufficiently 
large to reverse the sign of i>\ this is observed for some 
MSPs accelerated in the gravitational potential of globular 
clusters. 
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l2005l : lBovle fc Buonannoll2007h . For these sources, the 
strain spectrum is expected to be somewhat steeper 
with a ~ — 1. 

Observed pulsar periods are affected by many other 
factors. Pulsars in binary systems of course have their 
apparent period modulated by their binary motion and 
by relativistic effects. Intrinsic pulsar periods are not 
perfectly stable, with two main types of irregularity be- 
ing observed: apparently random period "noise" and 
sudden jumps in the period known as glitches. For- 
tunately, the s trength of these irregulari t ies is corre- 
lated with \i>\ llArzoumanian et al.l[l994l : IWang et alJ 
l200d : ISh annon fc Corded 120101 ) and so they are weak 
or undetec table in MSPs as th ese have very low val- 
ues of 1 1> | IjHobbs et alJl2010bh . Pulsar signals suffer 
frequency-dependent propagation delays in the inter- 
stellar medium, mainly due to dispersion and usually 
parameterised by the dispersion measure (DM), com- 
puted using 



where K = 2.410 x 10" 16 Hz -2 cm -3 pc s _1 and ti, t 2 
are pulse ToAs measured at radio frequencies fi , f 2 re- 
spectively. Because of the changing path to the pulsar 
resulting from the motion of the pulsar and the Earth 
in the Galaxy, these delay s are time- variabl e and must 
be continually monitored (|You et a l. 2007a) . Errors in 
the reference atomic timescale and in the ephemeris 
used for correction to the solar-system barycentre will 
directly affect the observed timing residuals. Finally, 
the data are affected by a variety of instrumental ef- 
fects such as receiver non-linearities, calibration errors 
and radio-frequency interference. 

These often unknown contributions to the timing 
residuals make it effectively impossible to detect GWs 
with observations of just a few pulsars. However, data 
from just one or two pulsars can be use d to place a 
limit o n the GW strength. As an example. iKaspi et al.l 
(|!994f ) used Arecibo observations of PSR B 1855+09 
and PSR B1937+21 over an eight-year data span to 
place an upper limit of about 6 x 10~ 8 on the en- 
ergy density of a stochastic GW background relative to 
the closure density of the Universe f2 gw at frequencies 
around 5 nHzQ 

For a positive detection of a stochastic background 
of GWs, its effects on the observed timing residuals 
must be separated from other contributions to those 
residuals. Fortunately, a method exists to achieve this 
separation. This method depends on making quasi- 
simultaneous observations of a number of pulsars that 
are distributed across the celestial sphere, thereby form- 
ing a "Pulsar Timing Array" (PTA). 

For an isotropic stochastic GW background, the 
expected correlation between residuals for pairs of pul- 
sars in different directions is a function only of the an - 
gle between the two pulsars (Hellings & Downs 1983). 
The precise dependence of the correlation on the sep- 
aration angl e depends on t he properties of the GWs 
themselves (|Lee et a l. 2008). According to standard 



3 The relative GW energy density is given by O gw (/) = 
ItFt/ 2 M/) 2 dPhinnevll20ul1;IJenet et al.|[200l) 



GR, GWs are spin-2 transverse-traceless waves. The 
correlation curve for such waves has a maximum of 
0.5 for pulsars which are close together on the sky 
(reduced from 1.0 because the same GW background 
passing over the pulsars produces a statistically equal 
but uncorrelated modulation in their residuals), goes 
negative for pulsars separated by about 90° and pos- 
itive again for 180° separation (see, e.g., Hobbs et 
al. 2009a). This well-defined "quadrupolar" signa- 
ture may be compared with the dipolar signature re- 
sulting from an error in the solar-system ephemeris, 
which effectively is an error in the assumed position of 
the Earth relative to the solar-system barycentre, or 
the monopolar signature resulting from time-standard 
errors which affect all pulsars equally. Other errors, 
for example, in the assumed parameters of a binary 
system, can be separated by their dependence on bi- 
nary orbital phase. Errors in interstellar corrections 
can be separated by their dependence on observing 
frequency. Careful attention to other corrections, for 
example, timescale transform ations, can ensure th at 
their uncertainty is negligible (Edwards et al]l200d\ 

The idea of observing a set of widely distributed 
pulsars, both to reduce statistical uncertainties and to 
enable the separation of timing pe rturbations as de- 
scribe d above, was introduced by iHellings fc Downs! 

(| 19831 ) and further discussed by iRoman i (1989) ( who 

coined the term "Pulsar Timing Array" ) and lFoster fc Backer! 
(|1990| ). 

The Parkes Pulsar Timing Array (PPTA) is an im- 
plementation of the PTA concept. It builds on ear- 
lier MSP ti ming efforts at the Parkes 64-m radio tele- 
scope (e.g.. iJohnston et al"1 [l993: Tosc ano et al.lll99^ : 
Ivan Straten et alj|200lft . The PPTA project is a col- 
laborative effort, primarily between groups at the Com- 
monwealth Scientific and Industrial Research Organ- 
isation (CSIRO) Astronomy and Space Science Divi- 
sion (CASS) and the Swinburne University of Technol- 
ogy with important contributions from astronomers at 
other institutions, both within Australia and interna- 
tionally. It uses the Parkes 64-m radio telescope to 
observe a sample of 20 MSPs at intervals of two to 
three weeks at three radio bands: 50 cm (centre ra- 
dio frequency f c ~ 700 MHz), 20 cm (~ 1400 MHz) 
and 10 cm (~ 3100 MHz). As well as the observa- 
tional program, the PPTA project has supported work 
on the analysis of pulsar timing data and its inter- 
pretation. A major component of this was the devel- 
opment of the T empo2 pulsar timin g analysis pack- 
age, d escri bed b y Hobbs et al] (|2006h . lEdwards et all 
( 2006) and iHobbs et all l|2008n . The PP TA has also 
suppo rted develop ment of the Psrchive plotan et ail 
2004b) and DSPSR (|van Straten fc Bailesll201ll ) pulsar 
data analysis packages. 

Detection of the GW backg round by a PTA was 
discussed by Ijenet etlil ] [|2005h . They showed that 
precision timing of at least 20 pulsars with ToA er- 
rors ~ 100 ns and a five-year data span was neces- 
sary for a positive detection of the expected stochas- 
tic GW background from binary super-massive black 
holes in distant galaxies. This work formed the ba- 
sis for the d e sign o f the PPTA observational program. 
iJenet et"al] (|2006r ) developed a method of analysing 
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PTA data sets to place a limit on the energy den- 
sity fi gw of the stochastic background based on the as- 
sumption of "white" or u ncorrelated t i ming residuals. 
They applied this to the iKaspi et al.l l| 19941) Arecibo 
data for PSR B1855+09 and Parkes observations of 
seven PPTA pulsars from early 2003 to mid-2006 to 



derive a limit on 
1. 9 x 

bv lKaspi et al 



n gw at a frequency of 1/(8 yr) of 



■y 

19941 1. lYardlev et all (|2010h discussed 
the sensitivity of a PTA to an isolated GW source, 
deriving the first realistic sensitivity curve for a PTA 
system. Using the method developed by IWen et all 
(120111 ) . lYardlev et al.l l|2010] ) also placed a sky -averaged 
constraint on the merger rate of black-hole binary sys- 
tems with a chirp mas a of io 10 m in galaxies with 

redshift z < 2 of less than one every 5 yr (see also 
Yardlev 2011). A method of detecting the stochastic 
GW background based on the Hellings and Downs cor- 
relation and considering the effects of irregular time 
sampling and spectral le akage on the detectio n sen- 
sitivity was developed by lYardlev et all (120111). This 
was ap plied to the Parkes o b serva tions of lVerbiest et all 
(|2008h and IVerbiest et"aU (|2009D showing that these 
data were consistent with a null result to 76% confi- 
dence. 

PTA observations of sufficient duration can in prin- 
ciple detect currently unknown solar-system objects, 
for example trans-Neptunian objects (TNOs). To date 
there has been no attempt to directly solve for an gen- 
eral dipole signature since the unknown dipole axis 
makes a general solution difficult. However it is possi- 
ble to test for predictable signatures such as errors in 
the m ass of the solar-system planets. I Champion et al.l 
(2010) used this method with the additional assump- 
tion that a change in the mass of a planet simply 
shifted the barycentre proportionally along the barycentre- 
planet axis, thereby introducing a sinusoidal varia- 
tion at the planetary orbital period into the observed 
residuals. They showed that the derived mass of the 
Jupiter sy stem was consistent with the value assumed 
in DE421 (jFolkner et al.ll2 008l to within 2 x 10" 10 M , 
a higher precision than that obtained from observa- 
tions of the Pioneer and Voyager spacecraft and con- 
sistent with the value obtained from Galileo. 

All ToAs are initially measured with reference to 
a local observatory clock. They are transferred to 
an international standard of time, for example, In- 
ternational Atomic Time TT(TAI), or the annual up- 
date of this, TT(BIPMxx), using clock offsets pub- 
lished by the BIPM0 Many authors have discussed 
the idea of establishing a "pulsar timescale" b ased on 
the rotation or orbital mot i on of pulsars (e.g ., iTavlorl 



1991 



2008) 



iPetit fc Tavellalll996l ; lllvasov et a"l]|1998l ; iRodinl 
A pulsar timescale is fundamentally different 
to timescales based on atomic frequency standards in 
that it is based on different physics - rotation of mas- 
sive bodies - and is largely independent of Earth-based 
effects, for example, seasonal variations. Another im- 



4 The chirp mass M c = (M 1 M 2 ) 3/5 (M 1 + M 2 )~ 1/5 , 
where Mi and M 2 are the masses of the two components of 
a binary system, is the mass parameter relevant to em ission 
of GWs from binary systems (see, e.g., Thornc 1987) 

5 See |http: / / www. bipm.org 



portant point is that MSPs will continue to spin for 
billions of years whereas individual atomic clocks have 
lifetimes measured in years or decades at best. Pul- 
sars may therefore provide a uniquely stable long-term 
standard of time. A pulsar timescale is essentially 
independent of the reference atomic timescale but it 
is not absolute and cannot measure linear frequency 
drifts of the reference atomic timescale. However, vari- 
ations of higher order can in principle be detected by 
comparison with a pulsar timescale. Two rec ent ef- 
forts at est ablishing a p u lsar t im escale a r e bv IRodinl 
(|201ll ) and iHobbs et all (|2012h . IRodinl (|201ll ) used 
seven years of archival timing data for six pulsars from 
the Kalyazin Radio Astronomy Observatory and an 
analysis method based on Wiener filtering to establish 
a pulsar times cale with stab i lity ~ 5 x 10 -14 over the 7- 
year interval. IHobbs et al. (2012) used Parkes timing 
data for 19 pulsars spanning up to 17 years, includ- 
ing nearly six years of PPTA data, to establish the 
pulsar timescale TT(PPTAll). This had a stability 
more than an order of magnitude better than that of 
IRodinl l|201ll ) and showed significant departures from 
TT(TAI). These deviations closely matched the dif- 
ferences between TT(BIPMll) and TT(TAI) over the 
same time intervals, thus demonstrating both that pul- 
sar timescales can be of comparable precision to the 
best atomic timescales and that the post-corrections 
used to form TT(BIPMll) do improve the stability of 
the timescale. 

PTA data sets have many other applications in- 
cluding detailed studies of the indi vidual pulsars and 
studies of the interstellar medium. lYou et all (|2007al ) 
used PPTA data sets t o study variat ions in inte r stellar 
dispersion. lYou et all (|2007bl) and lYou et all (|2012l ) 
studied the electron density and magnetic field in the 
solar wind using observations of several PPTA pulsars 
that have low ecliptic latitudes. Detailed studies of 
the 20cm polarisation and mean p ulse profile s of the 
PP TA pulsars were p resented by lYan et aD (|2011al ) 
and lYan et all (|2011bl ) investigated rotation measure 
variations for the PPTA pulsars, both short-term vari- 
ations due to changes in the E arth's ionosphere and 
long-term interstellar variations. Oslowski e t alj (|201ll ) 
investigated the ultimate limits to precision pulsar tim- 
ing in the case of high signal-to-noise ratio (S /N) ob- 
servations using J0437— 4715 as an example. 

This paper gives an overall description of the PPTA 
project which formally commenced in 2003 June . Pre- 
limina r y reports on the pr o ject were published bvlHobbs| 
' 2005 1 ) . Manchester! il2006ft iManchesterl (|2008l ). lHobbs et all 
2002) and lVerbiest et al.l (|2010T ) . n this paper, the re- 



ceiver instrumentation and signal processing systems 
used and (in some cases) developed for the PPTA project 
are described in 32 The observational strategy and 
details of the observed MSP sample are described in 
fj3]and the results obtained so far are described in ||4] 
Prospects for international collaborations with other 
PTA projects are discussed in iJS] along with future 
prospects for PTA projects, especially in the era of the 
Square Kilometre Array (SKA). The extended PPTA 
data set, formed by ad ding earlier Parkes timing d ata 
for the PPTA pulsars (jVerbiest et al.ll200Sl . 120091 ) to 
the PPTA data set is described in Appendix A. 
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2 Instrumentation and Signal 
Processing 

2.1 Telescope and receiver systems 

The PPTA project is based on observations made us- 
ing the Parkes 64-m radio telescope in New South 
Wales, Australia, with a number of receiver and back- 
end systems^ Parkes is at a latitude of —33° and so 
the Galactic Centre passes within a few degrees of the 
zenith. With its zenith-angle limit of about 59°, it has 
an effective northern declination limit of about +25°, 
and so can see more than two-thirds of the celestial 
sphere, including all of the rich inner Galaxy. As men- 
tioned in iJTJ observations are made in three different 
radio-frequency bands : 10cm, 20cm and 50cm. At 
20cm, observations are mostly made u sing the centre 
beam of the Parkes multibeam receiver l|Stavelev-Smith et 
1 19961 ') but occasionally with the "H-OH" receiver when 
the multibeam receiver is not available. Observations 
at 10cm and 50cm are ma de using the dual-ba nd coax- 
ial "10cm/50cm" receiver l|Granet et al.|[2005h . Details 
of these receivers are given in Table [1] where f c is the 
nominal band centre and S ay s is the system equivalent 
flux density. The original 50cm band suffered from the 
presence of digital television (DTV) signals transmit- 
ted from Mt Ulandra, about 200 km south of Parkes. 
By 2009, the number of DTV channels in the band had 
increased to five, covering more than half the band, 
and it was decided to move the band up in frequency 
to 700 - 770 MHz (labelled "40cm"). For convenience, 
the band label "50cm" implies both 40cm and 50cm 
data in the remainder of this paper and we refer to 
"three-band" data sets which in fact include data with 
four band designations. The availability of the various 
receiver packages is shown graphically in Figure [1] 

All systems receive orthogonal linear polarisations 
and (except for the 50cm receiver) have a calibration 
probe in the feed at 45° to the two signal probes. A 
pulsed broad-band noise source is used to inject a lin- 
early polarized calibration signal which is used to cal- 
ibrate the gain and differential phase of the two signal 
paths. Because of the coaxial nature of the 10/50cm 
feed, the 50cm receiver has four signal probes with op- 
posite signals being added in 180° hybrids to form the 
two orthogonal polarisation channels. The calibration 
signal is injected into directional couplers located be- 
tween the hybrids and the preamplifiers. Coupling be- 
tween the (nominally) orthogonal feeds is low (< —25 
db) for all receivers except the Multibeam receiver, 
where it is about —12 db. The spectrum of the in- 
jected calibration signals and the system equivalent 
flux density are calibrated in flux density units using 
observations of the strong radio source Hydra A (3C 
218) as is described further below f £|2.4[l . 

For the 20cm and 10cm systems, the signals are 
down-converted to baseband in the focus cabin. After 
transmission to the receiver control room, the signals 

6 We use the term "receiver" to mean all components of 
the receiving system from the feed to the digitiser inputs 
and "backend" to mean all components from the digitiser 
inputs to the data disks used for near-real-time storage of 
data. 



from all receivers are down-converted, band-limited, 

amplified and adjusted to the optimal signal l evel in a 

remo tely controllable down-conversion system (jGraves et al.l 
120001 ") for presentation to the backend digitiser systems. 

Observations are made under control of the Parkes 
Telescope Control System (TCS), a graphical interface 
which allows control of the telescope pointing and se- 
lection of the required receiver, backend system(s) and 
observation modes. PPTA observations are normally 
made under schedule control ensuring the correct se- 
lection of backend configuration and sequence of cali- 
bration and pulsar observations. 

2.2 On-line signal-processing systems 

A number of backend systems have been used in the 
course of the project. Their basic parameters are listed 
aUi Tabled The pulsar Wide Band Correlator (WBC) 
was a correlation spectrometer with an implementa- 
tion of the Fast Fourier Transform (FFT) algorithm 
on Canaris processor chips and 2-bit (3-level) digitiza- 
tion. The Parkes Digital FilterBank systems (PDFB1, 
PDFB2, PDFB3 and PDFB4) are (or were) implemen- 
tations of polyphase transforms using Field-Programmable 
Gate Array (FPGA) proc essors with 8-bit digitiz ation 
of the input signals (see iFerris fc Saunders] 120041 . for 
an example implementation and additional references). 
Polyphase filterbanks can be designed to have much 
superior channel bandpass characteristics compared to 
FFT spectrometers. For example, in our systems, the 
first sidelobes are —30 db down, compared to —12 db 
for an FFT-based system. Furthermore, the channel 
bandpass is much more rectangular. These properties 
lead to much improved isolation of narrow-band radio- 
frequency interference (RFI) and reduced signal loss. 
CPSR2 is a baseband recording system allowing coher- 
ent dedispersion of two 64- MHz- wide dual-polarisation 
data channels and APSR is its successor allowing co- 
herent dedispersion of dual-polarisation data at band- 
widths of up to 1 GHz. 

Regular observations with the WBC began in early 
2004, but because of various instrumental problems, 
high-quality data were not obtained until early 2005. 
The WBC was decommissioned in 2006 March. PDFB1 
was based on a commercial (Nallatech) signal-processing 
board and was in use from 2005 June to 2007 Decem- 
ber. PDFB2 was the first Parkes digital filterbank 
system to use processor boards developed at CASS, 
namely prototype board s for the Compact A rray Broad- 
band (CABB) system (| Wilson et al.l 1201 il l. PDFB2 
was commissioned with 1024-MHz bandwidth in June 
2007 and decommissioned in 2010 May. PDFB3 and 
PDFB4 use the final version of the CABB boards which 
have more powerful FPGA processors and have more 
on-board memory compared to those used for PDFB2. 
PDFB3 and PDFB4 were commissioned in 2008 July 
and remain in regular use. PDFB3 has two identical 
signal processing boards allowing simultaneous pro- 
cessing of up to four input signals and giving better 
performance in some configurations. Figure [5] shows a 
block diagram of the PDFB3 system. PDFB4 is simi- 
lar except that it has just one processor board and so 
cannot provide the inverse filter and subsequent pro- 
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Table 1: Receivers used for the PPTA project 



Receiver Band Freq. Range f c S sys 
(MHz) (MHz) (Jy) 



10/50cm 


50cm* 


650 - 


720 


685 


56 


10/50cm 


40cm 1 ' 


700 - 


770 


732 


62 


Multibeam 


20cm 


1230 - 


1530 


1369 


36 


H-OH 


20cm 


1200 - 


1800 


1433 


30 


10/50cm 


10cm 


2600 - 


3600 


3100 


50 



* Before 2009 July T After 2009 July 



cessing for the APSR system (described below). Fig- 
ure [1] shows the intervals when the various backend 
systems were in use for the PPTA project. 

The first five systems in Table [2] provide real-time 
computation of both the direct and cross products of 
the two complex baseband signals (A,B) giving full po- 
larisation information, i.e., A* A, B*B, Re(A*B) and 
Im(A*B). They also provide real-time folding of the 
channel and polarisation product data at the apparent 
pulsar period with the specified maximum number of 
bins across the pulse period. All systems are (or were 
in the case of the first three) controlled by a computer 
which obtains control data from TCS, including tele- 
scope, source and configuration data needed to sychro- 
nise the observations. The control computer operates 
on a basic cycle which is normally of 10 s duration. 
It provides displays of receiver bandpasses, digitiser 
levels and folded pulse profiles each cycle. 

The bandwidths given in Table[2]are maximum val- 
ues; all bandwidths less than the maximum by factors 
of two down to 8 MHz can be processed provided an 
appropriately band-limited signal is presented at the 
digitiser inputs. Observed bandwidths are normally 
64 MHz for the 50cm receiver, 256 MHz for the 20cm 
receivers and 1024 MHz for the 10cm receiver. The 
number of frequency channels and phase bins quoted in 
Table[2]are also maximum values in the sense that their 
product cannot exceed the product of the quoted val- 
ues; the number of channels and bins can be varied by 
factors of two within this limitation. Qualifications on 
this are that for the WBC, the quoted number of chan- 
nels is an absolute maximum and for PDFB1 the num- 
ber of channels could be changed downward only by 
factors of four. For PDFB2, PDFB3 and PDFB4, the 
maximum number of channels per polarisation product 
is 8192. In Table[5]the minimum folding period P m i n is 
given for the specified bandwidth, number of channels 
and number of bins. P m i n is (down to some limiting 
value) inversely proportional to the bandwidth and di- 
rectly proportional to the channel-bin product. For 
example, for a bandwidth of 256 MHz, 1024 bins and 
1024 channels, P min for PDFB3 and PDFB4 remains 
at 4.1 ms, but for 1024 MHz bandwidth, 512 bins and 
512 channels, it is just 0.26 ms. 

The WBC and PDFB1 systems used Tempo polyco 
files as a basis for predicting the topocentric folding 
period; PDFB2, PDFB3 and PDFB4 use(d ) Tempo2 
predi ction files which have greater precision (jHobbs et alJ 
120061 '). Pulsar astrometric parameters and barycen- 



tric period and binary data are obtained using PSR- 
CAT (Manchester et al. 2005) with a regularly updated 
database filefl An on-board "Pulsar Timing Unit" 
(PTU) increments the bin number at the rate required 
to maintain phase with the apparent pulse period, the 
latter being updated each cycle of the control com- 
puter, normally at 10-second intervals. The PTU and 
the input 8-bit digitisers use a 5 MHz signal locked to 
the Observatory H-maser as a primary reference; this is 
multiplied to 256 MHz and phase-locked to the Obser- 
vatory 1-second (1 PPS) signal for the PTU and then 
further multiplied to 2048 MHz for the digitisers that 
Nyquist-sample the 1024 MHz input bands. The PTU 
determines the lag between the zero of pulse phase at 
the start of an observation and the leading edge of the 
1 PPS signal with an uncertainty of 4 ns and sets the 
topocentric folding period or bin time each cycle with 
a precision typically better than 1 : 10 9 . 

Folded profiles for the four polarisation products 
for every frequency channel are transferred to the con- 
trol computer every cycle. They are integrated there 
for a "sub-integration" time, which is a multiple of the 
cycle time and normally 60 s, before being written to 
disk along with head er information in psrfits format 
(|Hotan et al.ll2004bl ')FI The psrfits output files also 
contain tables giving other information such as mean 
digitiser levels as a function of time, the distribution 
of digitiser counts, the receiver bandpass for each po- 
larisation channel, the pulsar parameters used for the 
prediction of the folding period and the predictor table 
used for the observation. 

CPSR2 (Caltech-Parkes-Swinburne Recorder, Ver- 
sion 2) was a baseband recording system which recorded 
two pairs of 64-MHz- wide baseband signals wit h 2-bit 
(4-level) digitization l|Bailesll2003 : lHotanll2007l '). Data 
were distributed to two primary processors and then 
to 28 secondary processors that performed coherent 
dedis persion of the baseband data ( Hankins fc Rickettl 
I1975T ). followed by folding at the apparent pulsar pe- 
riod to form 128-channel pulse profiles, typically with 
1024 pulse-phase bins per polarisation product. CPSR2 
was commissioned in August 2002 and decommissioned 
in July 2011. CPSR2 recorded data files every 16 s 
in "timer" format. These files were visually checked 
for RFI and then, if clean, summed to form 64-s sub- 
integrations. The sub-integration files for each obser- 
vation were then combined with additional header in- 



See http:/ /www. atnf.csiro.au/research/pulsar/psrcat 
See also http: / /www. atnf.csiro.au/research/pulsar/psrfits 
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formation to form a psrfits file. 

APSR (ATNF-Parkes-Swinburne Recorder) is a base- 
band recording system which uses PDFB3 for digiti- 
sation and signal conditioning (Figure [2}. It provides 
16 pairs of baseband signals covering a maximum to- 
tal bandwidth of 1024 MHz. Because of data transfer 
limitations, for 1024 MHz and 512 MHz total band- 
width, samples are truncated to two and four bits 
respectively; for smaller bandwidths, 8-bit data are 
recorded. Data are transferred via four 10 Gbit ether- 
net lines and a fast switch to a cluster of 16 dual-quad- 
core processors where coherent dedispersion and fold- 
ing a re carried out using DSPSR (jvan Straten fc Bailesl 
2011). The data are stored on disk at Parkes and then 
converted to psrfits format for subsequent process- 
ing. APSR has a web-based real-time monitoring and 
control system which communicates with TCS to syn- 
chronise the data recording. APSR commenced regu- 
lar operation in 2009. 

PDFB3 has an option for real-time rejection of 
radio-frequency interference (RFI) which can be used 
for both fold-mode and APSR observations. Two types 
of rejection are provided: a) time-domain clipping of 
impulsive broad-band interference and b) frequ ency- 
domain adaptive filtering o f quasi-stationary RFI l|Kesteven 
120051 ; iKesteven et alJl201Cf ). The latter requires a ref- 
erence signal containing the RFI which is fed to one of 
the second pair of digitiser inputs. Provided the RFI 
in the reference signal has sufficient S/N, the adap- 
tive filter removes the RFI from the two polarisation 
channels without disturbing the underlying astronomy 
signal. Its main application was to the original 50cm 
band (650 - 720 MHz) where digital television signals 
from transmitters on Mount Ulandra, located approxi- 
mately 200 km south of Parkes, were at significant lev- 
els. The reference signal was obtained from a 4-m re- 
flector with a vertically polarised dipole feed (to match 
the transmitted polarisation) directed at the horizon 
in the direction of Mount Ulandra. The same refer- 
ence signal cleans both polarisations of the astronomy 
signal, preserving the pulse polarisation. 

Quasi-real-time monitoring of pulsar profiles as a 
function of time, frequency and pulse phase together 
with input bandpass profiles is provided for all operat- 
ing backend systems by a web-based monitoring sys- 
temQ 

2.3 Calibration 

Calibration of the data is important to reduce system- 
atic errors associated with different bandpass gains and 
phases for the two polarisation channels, to place the 
Stokes parameters in a celestial reference frame, and 
to correct for the e ffects of cross-coupling in the feed 
van Stratenll200l. | 2006). Parameters describing the 
orientation of the signal and calibration probes rela- 
tive to the telescope a xis are stored in the mai n header 
of each psrfits file l|van Straten et alJl201Ch . Short 
(typically 2 min) observations of the pulsed calibration 
signal preceding (and sometimes following) each pul- 
sar observation are used to determine the instrumental 

9 http: / /pulseatparkes. atnf.csiro.au/dev/ 



gain and phase. The calibration data are then ap plied 
to the pulsar observations using the PsrchiveFi pro- 
gram PAC to flatten the bandpass and transform the 
polarisation products to Stokes parameters. Parallac- 
tic rotation is also corrected to place the Stokes pa- 
rameters in the celestial reference frame. 

For all systems, input signal levels are automati- 
cally adjusted to the optimal operating point (within 
±0.5 db) as part of the calibration procedure. For 
the PDFB systems, the operating point was chosen 
to give an rms variation of 10 digitiser counts which 
ensured linearity while preserving adequate headroom 
for strong signals. 

The pulsar data are placed on a flux density scale 
utilising observations of Hydra A, assumed to have a 
flux density of 43.1 Jy at 1400 MHz and a spectral in- 
dex of —0.91 over the PPTA frequency range. These 
observations are normally made once per session for all 
three bands and consist of a sequence of five 2-min cal- 
ibration observations at positions off north - on source 
- off south - on source - off north, where the off-source 
positions are 1° from the source position. The data 
are processed using the Psrchive program fluxcal 
to pr oduce flux calibration files for each band. These 
atrealsubsequently used by PAC to calibrate the pulse 
profiles in flux density units. 

A sequence of observations of PSR J0437— 4715 are 
made for each receiver system several times a year to 
calibrate the feed cross-coupling. Usually 8-10 ob- 
servations, each of 16-min duration, are made across 
the 10.5 hours that the pulsar is above the telescope 
horizon. These data are processed using the Psrchive 
program PCM with the "Measurement Eq uation Mod- 
elling" (MEM) method (|van Stratenl2004h to form "pcm" 
files. These can be used by PAC to correct the pulsar 
data files for the effects of feed cross-coupling. For 
this paper, cross-coupling corrections are only applied 
to data obtained using the 20cm multibeam receiver. 

2.4 Off-line Signal Processing 

All manipulation, visualisation and analysis of pulse- 
profile data is done using the Psrchive pulsar signal 
processing system. Data may be calibrated, shifted 
in pulse phase, summed in time, corrected for dis- 
persion delays and summed over frequency channels. 
Psrchive routines are also used for RFI excision and 
ToA estimation. A wide variety of different display for- 
mats is available, displaying the data as functions of 
pulse phase, time, frequency, polarisation parameters, 
etc. Programs to list file header and profile data in 
simple ascii formats are also provided. Figure [3] shows 
typical displays for PSR J1713+0747 from a standard 
64-min observation at 20cm with PDFB3. 

Pulse ToAs are important for many research areas 
including, of course, the PPTA project. While there 
are several options within Psrchive, ToAs and their 
uncertainties are obtaine d by perfor ming a Fourier- 
domain cross-correlation (|Tavlorlll99fjh of the observed 
pulse profile with a standard template for each pulsar. 
Noise-free standard templates were formed by inter- 

10 See http://psrchive.sourceforge.net 
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actively fitting scaled von Mises functions (using the 
Psrchive program paas) to a high S/N observed pro- 
file formed by adding many individual observations 
(cf., lYan et al.l l2011al >. The number of fitted func- 
tions varied depending on the complexity of the pro- 
file and was increased until the peak residual was less 
than or about three times the baseline rms deviation. 
For all pulsars except PSR J0437-4715, the profiles 
and template are total intensity or Stokes /; for PSR 
J0437— 4715 it is advantageous to fit to the invariant- 
interval profile, i.e., the profile of (I 2 — P 2 ) 1//2 , where I 
is the total-intensity Stokes parameter and P is the po- 
larised part of the signal, P = (Q 2 + U 2 + V 2 ) 1/2 , where 
Q, U and V are th e Stokes param eters describing the 
signal polarisation (Britton 2000). This largely avoids 
possible systematic errors associated with polarisation 
calibration of the data. Figure [4] shows mean observed 
pulse profiles, component von Mises functions and the 
template formed by summing these for three represen- 
tative PPTA pulsars at each of the three observing 
bands. For PSR J0437-4715 there are 13, 14 and 
17 von Mises components at 10cm, 20cm and 50cm 
respectively. Other pulsars have less than this, but 
never less than three (for J1909— 3744 at 10cm where 
the profile is narrow and has little structure) . The ToA 
reference phase is at the highest peak of the 20cm pro- 
file. Template profiles at 10cm and 50cm for a given 
pulsar were cross-correlated with the 20cm profile and 
aligned for maximum correlation. 

Each instrument has an effective signal delay (or 
advance) resulting from details of the time-tagging mech- 
anism and signal processing delays. These can be 
many tens of microseconds and must be calibrated be- 
fore data from different instruments are combined or 
compared. For the PDFB systems, delays were mea- 
sured by modulating the system noise with a PIN- 
diode in the signal line just before it enters the down- 
conversion system. The modulation signal was gener- 
ated by a programmable pulse generator and consisted 
of a pulse train (usually 6 or 12 pulses) with 40% duty 
cycle and total duration of about 0.95 ms, the first of 
which was triggered by the leading edge of the 1-ms 
pulse from the Observatory clock. This trigger is also 
synchronous with the leading edge of the 1-sec clock 
pulse which is used for the PDFB time tagging. This 
pulse train was observed with each instrument with a 
folding period of exactly n ms where n was typically 
2-5. Observation times were typically 2 min; several 
observations were made for each configuration, often 
in different sessions. 

A simulated pulse train matching the modulation 
signal was then convolved with the impulse response 
corresponding to the channel width of the particular 
PDFB configuration to produce a reference template 
for use with the Psrchive program pat used to pro- 
duce ToAs. An example of such a convolved template 
is given in Figure [5] The offset of the observed pulse 
profile relative to the template profile, modulo 1 ms, 
given by PAT is the instrumental delay for the partic- 
ular instrument and configuration. This measurement 
typically had an uncertainty of a few ns for a given ob- 
servation, but additional variable delays of 10 - 200 ns 
were revealed for some instruments. These variable de- 



lays, which occur apparently randomly and remain to 
be identified, may contribute a small amount of addi- 
tional effectively white noise to observed ToAs. Addi- 
tional constant offsets of 609 ns (estimated cable delay 
from the focus cabin to the point where the PIN mod- 
ulator was inserted) , 400 ns (instrumental delay in the 
pulse generator) and 30 ns (effective propagation delay 
from the focus cabin to the intersection of axes of the 
telescope) were subtracted from the measured delays 
to refer derived pulse ToAs to the intersection of the 
azimuth and elevation axes (the topocentric reference 
point) of the telescope. 

Correction for these instrumental delays allowed 
most of our observations to be placed on a common 
timescale. However, for some of the early instruments 
and configurations, these measurements were not made. 
In these cases, differential offsets from instruments with 
measured delays were determined by comparison of si- 
multaneous or contemporaneous ToAs for several of 
the more precisely timed pulsars. For a few of the 
most precisely timed pulsars, simultaneous measur- 
ments with different instruments revealed systematic 
ToA offsets, generally < 100 ns, which probably re- 
sult from a frequency- dependence of the relevant pulse 
profile. These were compensated for by a small phase 
rotation of the appropriate profile template. Measured 
delays are contained in a file which is automatically ac- 
cessed by the PPTA processing pipeline, the operation 
of which we discuss in £12.4. II Occasionally, file header 
parameters were incorrect, usually during the commis- 
sioning phase for a new instrument. These header er- 
rors are also corrected in the processing pipeline. 




0.2 0.4 0.6 0.8 

Pulse Phase 



Figure 5: Template profile used for measurement 
of instrumental delays for the PDFB3/4 configura- 
tion with 256 MHz bandwidth, 1024 channels and 
1024 profile bins. The rectangular input waveform 
has a leading edge at phase 0.0; the convolution 
for the finite channel bandwidth smears the edge 
transitions by a small amount. 

Observed ToAs are intially measured with refer- 
ence to the Parkes Observatory hydrogen maser fre- 
quency standard, UTC(PKS). These are referred to 
one of the international timescales published by the 
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Bureau International de Poids et Mesures (BIPM), for 
example, TT(TAI) or one of its retroactive revisions 
TT(BIPMxx), where 20xx is the year up to which th e 
timescale is computed (see, for example, Petit 2005). 
In this paper, we use TT(BIPMll), available at the 
BIPM ftp websiteE] Two methods of time transfer are 
currently available. The first is based on a Tac32Plus 
Global Positioning System (GPS) clock which directly 
gives UTC(GPS) - UTC(PKS) at 5-minute intervals. 
The BIPM publishes tables (in Circular T) from which 
daily values of TT(TAI) - UTC(GPS) and hence TT(TAI) 
- UTC(PKS) may be computed. The second system 
uses a GPS common-view link to UTC(AUS), operated 
by the National Measurement Institute in Sydney, giv- 
ing UTC(AUS) - UTC(PKS) also at 5-min intervals. 
Circular T also gives daily offsets of UTC - UTC(AUS) 
from which TT(TAI) - UTC(PKS) may be derived. 
Daily averages of these clock corrections are available 
to Tempo2. Typically the two modes of time transfer 
differ by a few nanoseconds after removal of a constant 
offset. 

Observed ToAs must be referred to the solar-system 
barycentre (SSB), assumed to be an inertial (unaccel- 
erated) reference frame, before comparison with pre- 
dicted arrival times based on a model for the pul- 
sar. This transformation uses the Jet Pr opulsion Lab- 
orato ry DE421 solar-system ephemeris (|Folkner et al.l 
2008) to correct for the motion of the Earth relative 
to the SSB, as well as oth er terms discussed in de- 
tail by Edward s et al.l (|2006f > . Variations in interstellar 
and solar System dispersion also affect observed ToAs 
and must be corrected for using measured or estimated 
dispersi on measures to re fer the ToAs to infinite fre- 
quency (|You et al.ll2007al lbTl. 

The PPTA project involves frequent observations 
over a long data span, generating large amounts of 
data, approximately two terabytes per year. Imme- 
diately after completion of an observation, data files 
are automatically transferred to archive d isks at CASS 
Headq uarters in Sydney. As described bv lHobbs et al.l 
(|201lT l , raw data files for most of these observations are 
publically available on the CSIRO Data Access Por- 
talQ which is part of the Australian National Data 
service after an embargo period of 18 months from 
the time of observation. 

2.4.1 Processing Pipeline 

As data files are transferred to Sydney, details of the 
observation including file name, pulsar name (or cali- 
bration identifier), telescope name, pointing right as- 
cension and declination, receiver, observing frequency, 
backend instrument and configuration, bandwidth, num- 
ber of channels, number of subintegrations, number of 
polarisations, project code and psrfits version num- 
ber are recorded in a mySQL table. Separate tables 
are used for calibration files, pulsar observation files, 
flux-calibration files, PCM-calibration files and profile 
template files. A processing script is then run to exe- 
cute the following sequence of commands for each new 

11 ftp://tai.bipm.org/TFG/TT(BIPM) 
12 https:/ /data. csiro.au/dap/ 
13 http:/ /www. ands.org.au 



calibration or observation file. Relevant Psrchive pro- 
grams are given in parentheses. 

• For pulsar files, the DM and pulsar parameter 
table are updated, (pam) 

• The observing band (10cm, 20cm or 50cm) is 
determined from the header data. Frequency 
ranges known to be contaminated by RFI for a 
given band and observation time are given zero 
weight. Band edges (5% of the bandwidth) are 
also given zero weight, mainly to remove aliased 
out-of-band signals. Pulse profile phase bins af- 
fected by undispersed transient RFI are set to a 
local mean, (paz) 

• Calibration files are checked for a bad first sub- 
integration (which occasionally occurs because 
of faulty synchronisation of the observation start) ; 
the sub-integration is given zero weight if faulty. 
(cleanCalFile) 

• For observation files, data are averaged in time 
to give files with eight sub-integrations. Cali- 
bration files are fully summed in time, (pam) 

• Header data known to be invalid are corrected. 
(fix_ppta_data) 

• Observation start times are adjusted to compen- 
sate for instrumental delays, thereby referring 
the start time to the telescope intersection of 
axes, (dlyfix) 

• Profiles are calibrated, correcting for instrumen- 
tal gain and phase and placed on a flux density 
scale, (pac) 

• This calibration step is repeated for 20cm data 
including the MEM calibration, (pac) 

• Profile data are summed to form 32 frequency 
channels and to form either the Stokes param- 
eters or (for PSR J0437-4715) the polarisation 
invariant interval. The resulting files are stored 
in the data archive and referenced in the database. 
(pam) 

• Profile data are fully summed across frequency 
and time and polarisations combined to form ei- 
ther Stokes I or invariant-interval profiles which 
are stored in the data archive and referenced in 
the database, (pam) 

• The fully summed profiles are cross-correlated 
with an appropriate template profile to form pulse 
ToAs which are stored in the database (pat) 

A set of scripts are available to obtain information 
from the database, for example, header data for a given 
observation or pulse ToAs for a given pulsar with or 
without the MEM calibration. The ToAs are provided 
in Tempo2 format and include flags for the observing 
band and receiver-backend system. Tempo2 can use 
these flags as well as observation times (in MJD), ob- 
servation frequencies or ToA uncertainties to select a 
particular subset of the observations using command- 
line arguments -PASS and -filter. More complex se- 
lection can be carried out within a select file that 
contains a list of filters. 
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Observations that are affected by remaining RFI, 
calibration or other instrumental errors are examined 
and corrected if possible. If uncorrectable, the affected 
observation is flagged as bad in the database and not 
included in subsequent processing or analysis. 

3 Observational Strategy 

Simulations of the predicted stochastic GW backg round 
from binary supermassive black holes in galaxies fl jaffc & 
2003l:IWvithe fc Loebll20o3 : lEnoki fc Nagashimall2007l : 
Sesana et al.ll2008l ) show that timing of about 20 pul- 
sars with weekly observations over a five-year data 
span with rms residuals of about 100 ns are needed for 
a sig nificant detection l|Jenet et al . 2005; Hobb s~et al.l 

2008) . For a background as described by Equation [TJ 
the (one-sided) power spectrum of the timing residuals 
is a power law given by 

A 2 f f \ 2a " 3 

(e.g.. iJenet et all 120061 ). For a = -2/3, the spec- 
tral exponent of the residual fluctuations is very steep, 
— 13/3. Consequently, the sensitivity of a PTA is greatly 
improved with longer data spans as the GW back- 
ground signal dominates the overall spectrum at low 
frequencies. As will be discussed more fully in iJS] be- 
low, the minimum detectable GW signal is roughly 
proportional to the amplitude of the timing residual 
fluctuation and therefore scales approximately as T" 13 / 6 ; 
where T is the data span. It also scales approximately 
as N , the number of pulsars in the PTA ( |Verbiest et al.l 

2009) . It is important to note that with significantly 
less than 20 pulsars, no matter how precise their ToAs, 
it would be impossible to make a significant detection 
of the expected stochastic GW background. This is 
primarily a consequence of the GW self-noise, that is, 
the ran dom noise introduced into the correlation sig- 
nature (jHellings fc Downdl 19831 : iHobbs et al.ll200Sft by 
the uncorrelated GWs passing over each pulsar. 

3.1 Sample Selection 

Only MSPs have the period stability and potential 
ToA precision to be useful for a PTA. About half of 
all known MSPs are located in globular clusters, but 
these pulsars are less useful partly because they are 
often very weak, but more importantly, because gravi- 
tational interactions with other cluster stars introduce 
additional perturbations to the observed pulsar pe- 
riod. Similarly, because of the relativistic perturba- 
tions, pulsars in very tight binary orbits are not ideal, 
especially if the orbit is eccentric. ToA uncertainties 
are approximately proportional to pulse width divided 
by the S/N. Consequently for timing array purposes, 
relatively strong Galactic-disk MSPs with short peri- 
ods and/or narrow pulses and either isolated pulsars 
or pulsars in wide binary orbits are preferred. Fig. [S] 
shows the distribution on the celestial sphere of MSPs 
suitable for timing arrays and those selected for the 



PPTA0 This figure illustrates the fact that at present 
there are few pulsars suitable for PTA observations 
that are inaccessible to the Parkes telescope. 

Table |3] lists the MSPs selected for the PPTA. 
Following the pulsar J2000 name, the pulsar period 
P, DM and orbital period Pb (if applicable) and the 
standard observation time are given. The next six 
columns give the mean and rms pulse flux densities 
(averaged over the observation time) for the 50cm (700 
MHz), 20cm (1400 MHz) and 10cm (3100 MHz) bands, 
Baffigf ectively. The final two columns give the mean 
pulse width at 50% of the peak level for the 20cm 
(1400 MHz) total intensity pulse profiles and its rms 
uncertainty. These flux densities and pulse widths 
are derived from PPTA observations as described in 
i)4.2l below. All of these pulsars, with one exception, 
are located in the Galactic disk. The exception, PSR 
J1824-2452A, is located in the globular cluster M28. 
It was included in the PPTA sample partly because it 
is relatively strong and can be accurately timed and 
partly to investigate the effects of period irregulari- 
ties on timing-array analyses. As well as the globular- 
cluster perturbations, this pulsar has relatively large 
DM variations. Furthermo re, a small period glitc h 
was reported for this pulsar (jCognard fc Backer 2004). 
As Fig. shows, the selected pulsars are widely dis- 
tributed on the celestial sphere and consequently pro- 
vide a good range of angular separations for the cor- 
relation analysis. 

There are a number of on-going searches for pulsars 
at various observatories around the world. When an 
MSP suitable for PTA observations is discovered and 
its basic parameters measured with sufficient accuracy, 
consideration is given to including it in the PTA obser- 
vations. From time to time, consideration is also given 
to dropping one or more of the less er-performing pul - 
sars. Specifically, PSR J2241-5236 (|Keith et alfcOllh 
was added to the PPTA obse rving schedule on 2010 
February 9, PSR J1017-7156 (|Keith et alj|2012af ) was 
added on 2010 September 7 and PSR J1732-5049 was 
effectively dropped from the schedule in 2011 January. 
Results for the two new pulsars are not discussed in 
this paper. 

The PPTA sample of MSPs has been regularly 
observed with good quality data at the three bands, 
50cm, 20cm and 10cm, since early 2005. In this pa- 
per we report on data obtained between 2005 March 1 
(MJD 53430) and 2011 February 28 (MJD 55620). Ob- 
servation sessions are typically 2-3 days in duration 
and at intervals of 2 - 3 weeks. For PSRs J1857+0943 
and J1939+2134, shorter observation times were cho- 
sen (Table [3} since these are northern sources which 
are being monitored at other observatories. PSR J2124— 
has a timescale for diffractive scintillation at 1400 MHz 
which is longer than the observation time and hence 
the pulsar is often not visible. If there is sufficient 
observing time, having two shorter observations sep- 
arated by a time long compared to the scintillation 
timescale increases the probability of obtaining a good 
ToA. Generally, low-DM pulsars tend to be more af- 
fected by scintillation since the scintillation bandwidth 
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is comparable to or larger than the observed band- 
width. Observers can terminate an observation pre- 
maturely if there is a low probability of obtaining a 
satisfactory ToA from the full observation time. This 
occurred occasionally, mostly for 20cm observations as 
the scintillation patterns are effectively uncorrelated 
for the 10cm and 50cm bands. 



3.2 Dispersion Corrections 

For the PPTA pulsars, DM variatio ns of up to 0.005 
cm -3 pc per year have been observed l|You et al.l2007ah . 
Since a DM variation of this size introduces a variable 
time delay of order 10 /is in a 20cm ToA it is clearly 
neccessary to correct for these variations. Timing anal- 
yses normally convert measured ToAs to infinite fre- 
quency using an estimate of the DM. If one is con- 
cerned about DM variations, this DM estimate may 
be obtained simultaneously with the timing analysis 
provided multi-frequency data sets are available. To 
illustrate the relevant issues we simplify the problem 
by assuming that the DM is measured using ToAs ti 
and t2 at two frequencies /i and /a, where /i is the 
primary frequency, i.e., the one with the best qual- 
ity ToAs. We use this DM to compute the infinite- 
frequency ToA ti jCX) corresponding to t\\ 



= ti(l+F)-t 2 F 



(4) 



where F = I/K/1//2) 2 - 1]. To ensure that the DM 
correction does not add significant noise to ti l00 , we 
require that the uncertainty in fe, 8t2, be less than 
[(1 + F)/F]Sti. For the PPTA observations, in the 
case where 20cm-band observations are DM-corrected 
using 50cm-band data, /1//2 ~ 2.0 and F ~ 1/3, and 
so the 50cm ToAs need to have an uncertainty less than 
four times that of the 20cm ToAs. For the 10cm/50cm 
combination, /1//2 ~ 4.4, F ^ 0.054 and the 50cm 
ToAs need to have an uncertainty less than 18 times 
that of the 10cm ToAs. The benefits of a large /1//2 
ratio are obvious. 

For the PPTA, the 50cm receiver is less sensitive 
than either of the 20cm or 10cm systems; the ratios 
of B 1 ^ 2 /S sys for 50cm:20cm:10cm are approximately 
(from Table Q] where B is the bandwidth) 1.0:3.3:4.8. 
The relatively lower sensitivity of the 50cm system is 
largely compensated for by the typically steep spec- 
tral index of pulsars. From Table the mean flux- 
density ratios across the 20 PPTA pulsars are {5700 : 
5i4oo) ~ 3.0 and (S700 : S3100) ~ 17.5, correspond- 
ing to mean spectral indices of —1.57 and —1.93 re- 
spectively. These ratios suggest that the 50cm ToA 
uncertainties should not contribute significantly to the 
uncertainty of the DM-corrected 20cm and 10cm ToAs. 
However, other factors also need to be considered. RFI 
is generally more of a problem at lower frequencies and 
so affects the 50cm data more than the other bands. 
Pulse widths are generally larger at lower frequencies, 
either because of an intrinsic frequency dependence of 
the pulse profiles or because of interstellar scattering. 
Some pulsars have flatter than average spectra and 
so the low-frequency ToAs have relatively greater un- 
certainties. Finally, ToA variations not described by 



the dispersion relation (Equation [2]) are sometimes ob- 
served at low frequencies. As a consequence of these 
factors, correction for DM variations is not always ad- 
va ntageous. These issu es are discussed in more detail 
bv lKeith et all (|2012bh . 



4 Results and Discussion 
4.1 Timing Data Sets 

ToAs produced by the processing pipeline described 
in M2.4.1I are stored in the mySQL database and are 
available for various timing analyses. As mentioned 
above, it is common for two or more backend instru- 
ments to simultaneously process the data from a given 
observation. ToAs based on different backend sys- 
tems processing the same input data are not indepen- 
dent and, for each observation, only the ToA with the 
smallest uncertainty is retained for subsequent analy- 
sis. For 50cm observations, the best quality data are 
always obtained with the coherently dedispersing in- 
struments, initially CPSR2 and later APSR. For the 
higher-frequency bands, the best data are normally 
obtained using the PDFB systems which have wider 
bandwidths. Exceptions are the high DM/P pulsars 
J1824-2452A and J1939+2134 at 20cm, for which the 
best results are usually obtained from APSR. 

For some pulsars the MEM calibration made lit- 
tle difference to the reduced \ 2 °f the timing solution 
whereas for others it resulted in a large improvement. 
For example, for PSR J1744-1134, the uncorrected 
20cm rms timing residual is 0.50 pis and the reduced 
X 2 is 11.0; for the corrected data the corresponding 
numbers are 0.32 u.s and 4.77. For simplicity, where 
the correction made little difference, the uncorrected 
data were used. 

Three of our pulsars have ecliptic latitudes less 
than 5°: PSRs J1022+1001 (-0.06°), J1730-2304 (0.19°) 
and J1824-2452A (-1.55°). ToAs from these pulsars 
are significantly delayed by the solar wind when the 
path to them is close to the Sun. Standard timing 
programs such as Tempo2 include a correction for this 
dispersive delay. However, observations show that the 
actual delay varies b y a factor of two or m ore from 
one year to the next (|You et al.ll2007bl . I2012T) and this 
variation is normally not modelled. The effect of this 
on our results is discussed in $4. 41 below. 

Figure [7] shows timing residuals for the final three- 
band data sets for four of the PPTA pulsars. These 
illustrate the relative timing precisions obtained in the 
different bands for different pulsars, mainly depending 
on the pulsar spectral index and the effect of DM vari- 
ations on t he timing PSR J 1045-4509 has large DM 
variations l|You et al.ll2007aT l and so there are signifi- 
cant residual variations that approximately follow the 
f~ 2 DM delay dependence. For PSR J0613-0200 on 
the other hand, the DM-related variations are small. It 
is also evident that there are systematic offsets between 
the ToAs in the different bands for a given pulsar. 
This is because the cross-correlation template align- 
ment procedure described in i)2.4l only approximates 
the actual frequency dependence of the profile. 
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Tables|4]-[6]summarise our timing data sets for the 
three bands. In each table, after the pulsar name, the 
next three columns are the full data span, the number 
of ToAs and the median ToA uncertainty, respectively. 
To give an indication of our current capabilities, the 
next four columns give results based on the last year of 
the PPTA data sets, that is, from 2010 March to 2011 
February inclusive. The columns are respectively, the 
number of ToAs in the one-year data set, the weighted 
mean timing residual and the corresponding reduced 
X 2 , and the median ToA uncertainty for the one- year 
set. Only the pulse frequency v and its first two deriva- 
tives, i> and v were fitted. The remaining parameters, 
excepting the DM offsets which were set to zero, were 
held at values obtained from the full three-band analy- 
sis described in 'A'i. 21 above. In a few cases, as noted in 
the tables, the data spans were varied to give parame- 
ters which more closely represented the system perfor- 
mance. For PSRs J1732-5049 and J1824-2452A, the 
data spans were increased to get a sufficient number 
of data points, and for PSR J1939+2134 they were de- 
creased to reduce the effect of the large DM variations 
on the rms residuals and reduced x 2 values. 

ToA uncertainties are as computed by the template 
fitting pr ogram, i n our case the Psrchive routine PAT 
using the Tavlog (|l990h algorithm. No scaling or bias- 
ing ("EFAC" or "EQUAD") factors have been applied. 
The tables show that in most cases the median ToA 
uncertainties are greater for the full six-year data sets 
than they are for the one-year data sets. This is es- 
pecially true for pulsars with high DM/P and mostly 
results from the improvement in back-end systems over 
the six years. 

It is notable that many of the reduced % 2 values 
for the one-year data sets are close to 1.0, especially 
for the 10cm data sets and for the weaker pulsars. 
This indicates that the computed ToA uncertainties 
are generally accurate. In some cases, the reduced % 2 
values are less than 1.0. This primarily occurs in pul- 
sars and bands that have large intensity fluctuations 
due to interstellar scintillation and hence a wide range 
of weights in the least-squares fit. This effectively re- 
duces the number of degrees of freedom for the fit and 
can result in an under-estimate of the corresponding 
rms residual. However, the minimum reduced x 2 value 
(0.58 for PSR J1024-0719 at 50cm) is still reasonably 
close to 1.0 and so the effect is not very significant. 

For the stronger pulsars, at 20cm and 50cm espe- 
cially, the reduced x 2 values tend to be larger than 
1.0, indicating short-term timing noise in excess of 
that expected from the ToA uncertainty. Possible rea- 
sons for this include residual RFI, especially at the 
lower frequencies, short-term variations in interstellar 
dispersion or scattering and short-term variations in 
the intrinsic profile shape (commonly known as "pulse 
jitter"). In some cases, e.g., PSRs J0437— 4715 and 
J1939+2134, the median ToA uncertainties are very 
small, < 30ns, and so reveal perturbations that are not 
obvious in the weaker pulsars. For PSR J0437— 4715 
it is likely t hat most of addition al scatter results from 
pulse jitter (jOsIowski et al.ll201lf h Errors in the forma- 
tion of the invariant interval used for timing this pul- 
sar are also possible. Although the invariant interval 



is no minally unaffected by calibration errors l|Brittonl 
2000), various second-order effects may be important. 
The invariant interval is essentially the difference be- 
tween Stokes I and the polarised component P. P is 
a positive-definite quantity and suffers a noise bias in 
the same way as the linearly polarized co mponent L — 
(Q 2 + U 2 ) 1/2 l|Lorimer fc Kramer! 120051 ) . This means 
that the invariant-interval profile is dependent on the 
S/N of the polarisation profiles used to form it. To 
minimise this effect, we summed the PSR J0437— 4715 
profiles to eight sub- integrations and 32 frequency chan- 
nels before forming the invariant intervals. Because of 
this issue, we did not use invariant-interval profiles for 
any of the other PPTA pulsars. 

PSR J1022+1001 is an interesting case. iKramer et alJ 
(1999) found variations in the relative amplitudes of 
profile pulse components on timescales of order hours 
and showed that such changes would induce offsets in 
derived ToAs. Their timing solutions had rms resid- 
uals of 20 - 100 /xs, substantially more than expected 
from random baseline noise, and they attrib uted the 
excess n oise to the profile variations. However. iHotan et al.l 
(|2004al 'l used carefully calibrated CPSR2 data recorded 
at Parkes and obtained an rms timing residual of just 
2.27 /xs in 5-minute integrations. They found no evi- 
dence for significant variations in profile shape on this 
or longer timescales. The 1-year 20cm MEM-calibrated 
PPTA data set has an rms timing residual of 0.66 ps 
(Table[5]). This shows that the effect of any profile vari- 
ations is very s mall, much less than those observed by 
IKramer et al.1 (|l999h . Never-the-less, the reduced x 2 
of 3.78 indicates that systematic ToA offsets are still 
present and that these have a timescale of the order 
of hours. As Figure H shows, PSR J1022+1001 has a 
huge frequency-dependent variation of profile shape. 
It also scintillates strongly, especially at 20cm (Ta- 
ble [3j). The pulse signal is therefore likely to be strong 
in different parts of the band at different times, re- 
sulting in a variable bias to the measured ToAs. This 
may dominate the ToA scatter that we observe. It 
is also possible that intrinsic temporal profile shape 
variations pla y a role, but a much smaller one than 
suggested by IKramer et al.l l|l999h . The pulse pro- 
file of PSR J1022+1001 has very high linear polari- 
sation, nearly 100% in the trailing half (jKramer et al.l 
119991 ; lYan et alJfeoild ). It seems most probable that 
the large profile an d timing variations observed by 
IKramer et al.l (|1999T ) resulted from errors in polarisa- 
tion calibration. 

Strong scintillation coupled with a frequency-dependent 
pulse profile may play a significant role for several of 
our pulsars. This effect could be largely eliminated 
by use of frequency-dependent templates, but it re- 
quires a careful study to assess its importance to the 
short-term timing noise relative to other contributions. 
For example, for PSRs J1824-2452A and J1939+2134 
there are large DM variations, the effects of which are 
not absorbed by the v—v—v fitting. Residual RFI also 
may make a significant contribution to the observed 
ToA scatter, especially at 50cm. Finally, instrumen- 
tal errors and deficiencies in signal-processing methods 
may contribute. These poorly understood contribu- 
tions to our timing residuals are the subject of contin- 
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uing investigation. With an improved understanding, 
we should be able to reduce their effect and hence ob- 
tain more accurate data sets. 

4.2 Flux Densities and Pulse Widths 

The mean and rms flux densities given in Table were 
computed by taking PDFB observations (which have 
a more reliable flux-density calibration) corresponding 
to ToAs from 2010 January to the end of the data set. 
The flux densities were averaged in time over each ob- 
servation and in frequency over each band. For PSR 
J0437— 4715, the total-intensity (Stokes I) profile was 
used; the invariant-interval mean flux densities are ap- 
proximately 60%, 85% and 80% of the Stokes I val- 
ues at 10cm, 20cm and 50cm, respectively. The frac- 
tional rms variations (or modulation indices) of the 
measured flux densities for each band are often large, 
especially for the low-DM pulsars. As mentioned in 
f° r such pulsars, observations when the pulsar had 
low flux density were sometimes terminated early, es- 
pecially at 20cm. The 20cm flux densities quoted in 
Tableware consequently somewhat biased to high val- 
ues, but this effect is small, typically a few percent or 
less. 

Table [3] also gives mean pulse widths at 50% of the 
profile peak and the estimated uncertainty in these 
widths. These results are averages of widths measured 
from the observation-averaged 20cm total-intensity pro- 
files. For the PSR J0437— 4715 invariant-interval pro- 
files, the mean width is 0.1370 ± 0.0004 ms, approxi- 
mately 3% smaller than the Stokes I width. 

4.3 Dispersion Variations 

DM variations have been determined for all 20 PPTA 
pu lsars using the three -band data as described in detail 
bv lKeith et all (|2012bh . Here we give a brief summary 
of the process. Tempo2 allows fitting to multiband 
data sets for a set of offsets from the nominal DM at 
specified times. The offsets are constrained to have 
zero mean and the average DM is obtained from the 
same fit. The sample interval was chosen to minimise 
the noise contributed to the DM-corrected residuals by 
estimation errors in the DM offsets. The optimal in- 
terval for each pulsar was taken to be the inverse of the 
"corner" frequency at which the power in the DM vari- 
ations equals the white noise level in the "best" data 
set (see below). To obtain DM offsets at other times, 
Tempo2 interpolates between the sample values. Sim- 
ply fitting for the DM offsets in this way will absorb 
some of the frequency-independent variations that are 
crucial to PTA objectives. To overcome this prob- 
lem, we chose to additionally fit for a "common-mode" 
frequency-independent offset sampled at the same in- 
tervals as the DM offsets. This common-mode signal 
can be used for PTA purposes, but for this paper, we 
only use it to provide an unbiased estimate of the DM 
variations. 

Because it is impossible to align template profiles 
for the different bands with sufficient accuracy, it is 
also necessary to simultaneously fit for offsets between 
bands (with one chosen as reference). A single offset 



per band pair is fitted. In the fit, To As are weighted 
by the inverse square of their uncertainty. Inclusion 
of other free parameters in the timing model depends 
on ToA precisions and the physical properties of the 
pulsar, such as its distance an d the nature of th e binary 
system. Cholesky whitening ( Coles et al. 20l"lh is used 
to properly account for the correlation in the residuals 
and thus to provide reliable error estimates for the 
DM variations and other parameters. Figure [8] shows 
the DM variations derived from the ToAs shown in 
Figure [7] using the above procedure. 

4.4 Single-band Corrected Timing Data 
Sets 

As a basis for applications of the PPTA data set, ToAs 
for the band having the lowest over-all rms timing 
residual, either with or without the MEM calibration 
and with or without correction for DM variations, were 
selected. Table [7| gives this "best" band, the correction 
procedure adopted and the data span in years for each 
PPTA pulsar. Figure [9] shows the profile templates 
used for the best instrument of the best band. The 
best instrument was PDFB4 in all cases except two. 
For PSRs J1824-2452A and J1939+2134 APSR was 
used because their high DM/P requires the coherent 
dedispersion provided by this instrument. The tem- 
plate reference phase is shown on each plot; knowledge 
of this allows comparison of the absolute pulse ToAs 
from the PPTA with data from other telescopes. For 
the 3-band Cholesky fit described in N4.3I the pulsar 
parameters, including pulse frequency v and z>, inter- 
band jumps and DM corrections were simultaneously 
fitted. The fifth column of Table [7] gives the number 
of pulsar parameters fitted and the sixth column gives, 
where applicable, the averaging interval for the DM 
corrections. The seventh and eighth columns give the 
rms timing residual and reduced x 2 after a fit of just 
v and v to the best-band data, with all other param- 
eters, including DM offsets where applied, held fixed 
at the values resulting from the 3-band fit. Since the 
best-band data sets have just a single band, there are 
of course no interband jumps. Figure [10] shows the 
post-fit timing residuals for these best-band fits. 

For three of the 20 pulsars we obtain our best tim- 
ing by using lOcm-band data. There are two main 
reasons why higher frequencies will give more precise 
ToAs. The first is that profile components generally 
tend to be narrower at higher frequencies, leading to 
improved timing precision. The second and usually 
more important reason is that interstellar propaga- 
tion effects are less significant at higher frequencies. 
Dispersion delays vary as f~ 2 and scintillation effects 
more like f~ 4 . We attempt to correct for DM varia- 
tions, but these corrections are always imprecise and 
the error induced in the residuals follows the same f~ 2 
dependence. 

Scintillation effects, especially those due to refrac- 
tive scintillatio n, are much more difficult to measure 
and correct for (Hcmbcrger fc Stinebrin e 2008; Coles et all 
l2010l ; lDemorestl201ll '). This is especially true for MSPs 
and up to now such techniques have not been success- 
fully applied to this class of pulsar. Scattering delays 
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have long-term variations as the path to the pulsar 
moves through the interstellar medium because of the 
relative transverse motion of the scattering medium, 
the pulsar and the Earth. However, diffractive scin- 
tillation also introduces short-term variations. Mean 
pulse profiles for most pulsars are significantly fre- 
quency dependent. Consequently, changes in the rela- 
tive pulsar flux density in different parts of the band 
will directly result in variable ToA offsets. This will 
vary on the timescale for diffractive scintillation. As 
indicated by the large rms variations in observed flux 
density given in Tabled this is either comparable to or 
larger than the observation time for many of the PPTA 
pulsars and bands. Refractive scintillation also intro- 
duces a variable delay due to variations in the effective 
path length as different parts of the scattering screen 
are focussed on the Earth. This is a broad-band effect 
and is not accompanied by any significant change in 
profile shape. It is therefore very difficult to predict 
and hence to correct for. 

The pulsars for which we obtain our best timing 
at 10cm are among the strongest PPTA pulsars for 
this band (Table [3}. This illustrates the fact that the 
primary reason that the 20cm band is superior to the 
10cm band for most of our pulsars, despite the fac- 
tors discussed above, is simply S/N. Because of the 
relatively steep spectrum of most pulsars, the low S/N 
obtained at 10cm outweighs the benefits of the higher 
frequency. This is also illustrated by the large rms tim- 
ing residuals we obtain at 10cm (Table [4]) for many of 
our pulsars. Because of these factors, it is likely that 
the best PTA-related use of more sensitive telescopes 
will be to allow more pulsars to be observed at higher 
radio frequencies. 

Nine of the 20 pulsars benefitted from the MEM 
calibration to the 20cm data and 14 of the 20 benefit- 
ted from correction for DM variations. For four pul- 
sars neither correction gave a significant improvement 
in the rms residuals. In principle, the MEM calibra- 
tion should always improve the stability of the profiles 
since it removes the pa rallactic angle dep endence of 
the calibrated profiles l|van Stratenl 12004). In prac- 
tice, the correction made little difference to the more 
weakly polarised pulsars and in these cases it was not 
applied. For PSR J0437-4715, results with the MEM 
calibration were inferior to those obtained by using the 
invariant interval. This pulsar is a special case in that 
the polarisation flips between orth ogonal states righ t 
at the peak of the profile (see, e.g.. lYan et al.ll2011al ). 
The ToAs for this pulsar are therefore extremely sen- 
sitive to any errors in the calibration which alter the 
intensity ratio of the orthogonal states at the pulse 
peak. The invariant interval is relatively immune to 
calibration errors and in this case gives the best re- 
sults. 

In most cases, rms residuals were reduced by ap- 
plying the correction for DM variations. For PSRs 
J1603-7202 and J1643-1224, systematic frequency- 
dependent residual variations that did not follow the 
f~ 2 DM law were observed. For both pulsars, applying 
the derived DM correction made the 20cm rms post-fit 
residuals significantly worse. In both cases, the inter- 
vals of non-DM residual variation were isolated and 



correlated in time over about six months, in 2006 for 
J1603-7202 and in 2010 for J1643-1224. It is possible 
that these were episodes of significant refractive delay 
but further investigation is needed. It is interesting 
to note that an extreme scattering event lasting three 
years and centred in 1998 was identified in flux-de nsity 
data for PSR J1643-1224 bv lMaitia etH] (|2003h . 

As mentioned in H4.1l above. solar- wind delays may 
be significant for at least three of the PPTA pulsars. 
Unmodelled delays are typically a few microseconds or 
less. However, since we do not always observe when a 
pulsar is close to the Sun, they affect just 5-10 ToAs 
even in the pulsars that pass within 5° of the Sun. 
We have chosen to include these ToAs in our data sets 
to allow further investigation of the solar-wind effects. 
The effect on the timing results is small. For the most- 
affected pulsar, PSR J1022+1001, eliminating ToAs 
for times when the pulsar was within 5° of the Sun 
from the single-band data set (Table 0) reduces the 
rms residual from 1.72 /xs to 1.64 [is and the reduced 
X 2 from 9.27 to 8.85. For all other pulsars, the effect 
is less significant. 

4.5 Long-term pulse frequency vari- 
ations 

The results presented in Figure [10] and comparison of 
the 6-year and 1-year rms timing residuals in Tables [1] 
-[7|show that, for a 6-year data span, about half of the 
PPTA pulsars have significant long-term fluctuations 
in their pulse frequency. These fluctuations may be 
quantified in different ways. A simple approach is to 
fit a pulse frequency second derivative (p) to each data 
set; this i s the essence of the timin g stability parame- 
ter A(t) (jArzoumanian et al.lll994T ) and quantifies the 
presence of a cubic term in the residuals. In Table [S] 
we present v values for the PPTA pulsars obtained in 
two different ways. The second column for each pulsar 
lists the v value obtained from a simultaneous fit of 
all parameters as described in £14.31 together with u, 
to the 3-band data sets. In the third column for each 
pulsar we list the v value obtained by adding it to the 
best-band fit as described in 34.41 that is, a simultane- 
ous fit of v, v and v to the best-band data set for each 
pulsar. 

In general, the v values from the 3-band fits have a 
lower significance than the values obtained by fitting to 
the best-band data alone. There are several reasons for 
this. Firstly, we are fitting for more parameters and 
some of these have significant covariances which are 
reflected in the derived uncertainties. Secondly, even 
though we try to minimise the effect by averaging over 
suitably chosen intervals, fitting and correcting for DM 
variations always adds effectively white noise to radio- 
frequency ToAs because the DM corrections are always 
uncertain at some level. Thirdly, our procedure for fit- 
ting v to the best-band data may under-estimate the 
true uncertainty in this parameter where DM correc- 
tions are applied since some correlation is introduced 
into the effective ToA errors because of the averaging 
of the DM corrections. This is not an issue for the 
seven pulsars where DM corrections are not applied to 
the best-band data sets (Table 0). In these cases the 
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fitting procedure is absolutely standard and parame- 
ters and their uncertainties should be reliable. 

In almost all cases, the v values derived by the two 
methods are in agreement within the combined un- 
certainties. A notable exception is PSR J0437— 4715. 
Because of its large flux density (Table [3} very pre- 
cise To As are in principle measurable for this pulsar. 
However, it is also true that observations of it reveal 
systematic errors most clearly. Unfortunately, as Fig- 
ure [10] shows, at its best band (10cm), data prior to 
mid-2006 were not of the same quality as subsequent 
data and therefore were not included in the data set. 
Earlier data at 20cm and 50cm are included in the 
3-band data set and so the two data sets are not di- 
rectly comparable, leading to the significantly different 
v values in Table [8] 

Half of the PPTA pulsars have a v value that is 
significant at the 3<r level based on the best-band fits. 
While there are some caveats about the uncertainty es- 
timations as discussed above, we believe the overall re- 
sult is solid. As an illustration, we take the important 
case of PSR J1909— 3744, where we find a significant v 
in contrast to other publish ed results l|Verbiest et al.l 
2009; [Demorest et all 120121 ') where no red noise was 
detected in data sets of similar duration[3 From Ta- 
ble [H] v values for the 3-band fit and the best-band fit 
are respectively 7.8 ±3.5 and 10.1 ± 2.1 x 10" 28 s" 3 . If 
we fit all 17 pulsar parameters plus v to the best-band 
data including the DM corrections held fixed, we ob- 
tain v = 14.6 ± 2.5 x 10" 28 s" 3 . If we omit the DM 
corrections, giving a standard multi-parameter fit to 
single-band data, we obtain v = 14.5±2.6 x 10 -28 s -3 . 
All of these values are consistent to 1.6cr or better and 
all the best-band fits give a value of v that is significant 
at about 5<j. 

These results represent the first detections of red 
timing noise in a significant sample of MSPs. Red tim- 
ing n oise has long been recognised in PSR J19 39+2134 



for which the ("-values are 0.51, 0.17, 1.68 and 0.54, 
respectively. Except for PSR J1824-2452A, the ob- 
served rms residual is within a factor of two of the 
prediction, which is satisfactory. There is some evi- 
dence that the model over-estimates the timing noise 
in PSR J1824-2452A Based on the spr e ad of tim- 
ing noise measured bv lShannon fc Cordesl (|2010l '). the 
value of £ for this pulsar is a la outlier. The model 
prediction for PSR J1824-2452A is high largely be- 
cause of the large This pulsar is unique in the 
PPTA sample in that it is associated with a globu- 
lar clus ter and furtherm ore is within the core of the 
cluster (|Lvne et al.lll987T l. suggesting that the large |z>| 
may partly result from acceleration of the pulsar in the 
cluster gravitational field. Consequently, no conclusion 
can yet be drawn about whether or not its timing noise 
is anomalously low. 

A more informative way to illustrate the pulse phase 
fluctuations for a given pulsar is to compute the power 
spectrum of the timing residuals. Figure [TTJ shows 
such spectra for the PPTA pulsars. These spectra 
were computed using a weighted least-squares fit to the 
timing residuals shown in F igure 1 101 after wh itening 
using the Cholesky method ijColes et al.ll201lf ). The 
line shown on each of the spectral plots has a slope 
of —13/3 for A g — 10 -15 (Equation [3]), a representa- 
tive value for the expected GW background from bi- 
nary superma ssive black holes in the cores of distant 
galaxies (e.g.. ISesana et all 12008'!. These plots illus- 
trate the wide differences in both red-noise and white- 
noise properties between our pulsars. Consistent with 
the residual plots of Figure[l0]and the rms residuals for 
the white noise given in Table [7] the white-noise power 
level differs by about four orders of magnitude from 
the "best" to the "worst" pulsars. More importantly, 
some pulsars have a much stronger red-noise compo- 
nent than others. PSR J1824-2452A stands out with 
the highest observed low-frequency noise level, most 



(agJKaspi et al.|1994|) and PSR J1824-2452 flHobbs et al.| pro bably as a result of the combined effects of intrinsic 
1200 " 



these pulsars are anomalous in having much 
stronger red noise than other MSPs. Indications of 
long-term correlated residua l fluctuations have been 
seen in MSPs J1713+0747 (ISplaver et al.ll2005l) and 
J0613-0200, J1024-0719, J1045-4509 ([Verbiest et alJ 
|2009|K but no v values were derived. 

[shannon fc Cordesl (|2010h chose to quantify the 
noise properties of pulsars using the rms timing residu- 
als rather than v. They modeled the rms timing resid- 
uals as a function of v, v and data span for different 
samples of pulsars. We have compared the rms tim- 
ing residuals in Table [7] with the values predicted by 
their "CP+MSP" model, which is based on a fit (in log 
space) to published rms timing residuals for both nor- 
mal (non-millisecond) pulsars and MSPs. They define 
a parameter £ which is the observed rms timing resid- 
ual after fitting for v and v divided by the predicted 
value. For the PPTA pulsars, most C, values are of or- 
der ten or greater, indicating that the observed rms 
residual is dominated by white noise. Exceptions are 
PSRs J0437-4715, J1824-2452A, J1909-3744 and J1939+2 



15 In fact, in his thesis Vcrbicst (2009), did quote a v for 
PSR J1909-3744, (11 ±4) X 10" 28 s -3 , which is consistent 
with the value in Table [8] 



period irregularities and the varying spatial accelera- 
tions due to gravitational interactions with other stars 
in the globular cluster. PSR J1939+2134 also stands 
out, but mainly because of its very precise To As; the 
absolute level of the red noise in this pulsar is compara- 
ble to that of several others, e.g., PSRs J1045-4509, 
J1603-7202 and J1643-1224. The ten MSPs which 
have v significant at 3<r or greater in the best-band fit 
(Table [5]) are marked with a * after the name in Fig- 
ure [TTJ Generally in these cases the power level in the 
lowest-frequency bins is above the white-noise level. 

It is important to note that, for the best pulsars, 
e.g., PSRs J0437-4715, J1713+0747, J1857+0943 and 
J1909— 3744, the observed power level at the lowest 
frequency (~ 0.17 yr" 1 ) is already less than the power 
expected from a GW background with A g — 10 -15 . 
However, there are some important caveats here. The 
GW background line is the average for an infinite num- 
ber of universes. Any given realisation (e.g., our own 
verse) may have a higher or lower spectrum. In 
ing a limit on the actual GW backgr ound, it is 
standard practice (e.g., iJenet et all 1200a ) to take a 
value for which 95% of the realisations have a detec- 
tion statistic that is greater than the statistic for the 
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observed spectrum. Also, comparison of the GW pre- 
diction with the measured spectra assumes that the 
Cholesky method compensates for the effects of fitting 
for v and v on the post-fit residuals and their spec- 
trum. The accuracy of this compensation depends on 
the accuracy of the estimation of the red-noise spectral 
model and this is difficult when the red and white noise 
are of comparable amplitude. Despite these caveats, 
we can safely conclude that the low-frequency noise 
level seen in our best pulsars is already close to the 
average level expected for an A s = fCP 15 stochastic 
background. 

Conversely, it is clear that the red noise in the 
"worst" of our pulsars cannot result from the GW 
background. If these signals were from a GW back- 
ground, they should be present in the timing residuals 
of all pulsars and they are not. In fact, these red sig- 
nals cannot result from any "global" process such as 
timescale or solar-system ephemeris errors which affect 
all pulsars. This is true even if there are geometric fac- 
tors that depend on the pulsar and/or source position, 
since it is very unlikely that our best pulsars would all 
happen to be located near the nulls of the geometric 
pattern. The most likely sources of pulsar-specific red 
noise are intrinsic timing irregularities and uncorrected 
interstellar variations. 

These results lead to a key question for PTA projects: 
to what extent will achievement of PTA goals be af- 
fected by intrinsic low-frequency pulsar period noise 
as data spans increase? Intrinsic fluctuations are of 
course uncorrelated between different pulsars and so, 
with a big enough pulsar sample, they will not pre- 
clude the reaching of PTA objectives which inherently 
depend on the detection of correlated signals. However 
there is no doubt that they make the task much more 
difficult, especially if the intrinsic noise is of greater 
amplitude than the target signal. 

The degree to which extension of data spans helps 
GW detection efforts depends critically on the spec- 
tral slope of the intrinsic red noise. As Table [8] and 
Figure [TTJ show, about half of the PPTA pulsars have 
detectable red timing noise. Estimation of the spectral 
properties of the red noise is difficult as, with a cou- 
ple of notable exceptions, viz, PSRs JI824— 2452A and 
JI939+2I34, the red signal is only marginally above 
the white noise level even at the lowest observed fre- 
quencies. Never-the-less, in order to make a first esti- 
mate of the spectral properties, we fitted a power law 
plus a constant (white) spectral model to the spectra 
shown in Figure If f I Only four pulsars (Jf 024— 0719, 
Jf 643-1224, JI824-2452A and Jf 939+2134) had power- 
law spectral indices greater than 2.5 and for only one of 
these (Jf 024— 0719) was the spectral index comparable 
to that of the expected GW spectrum. We emphasise 
that these estimates have considerable uncertainty, but 
it is clear that, for the majority of the PPTA MSPs, 
any red timing noise detected so far has a spectrum 
that is much flatter than the expected GW spectrum. 
We cannot preclude the presence of a currently invisi- 
ble red (non-GW) signal with a steep spectrum. How- 
ever, with that proviso and given current models for 
the origin of the GW background, we can assert that 
within the next decade a GW signal should be domi- 



nating the residuals of several of our pulsars and hence 
that a signficant detection of that GW signal can be 
expected. 

4.6 The Extended PPTA Data Sets 

The most important of the PTA objectives described 
in iJJJare for detection of phenomena that generate low- 
frequency fluctuations in residual time series, i.e., the 
spectra of the expected signals are very red. Conse- 
quently, it is highly advantageous to utilise data sets 
with long spans when searching for these phenomena; 
this is discussed more quantitatively in fJS]below for the 
case of GW detection. The PPTA data sets for the 20 
observed pulsars have data spans of approximately six 
years. However Parkes timing data were obtained for 
most of these pulsars for up to 11 years prior to th e 
start of the PPTA project (|Verbiest et al.ll2008l . [20091 '). 
In Appendix A, we present a re-analysis of these data 
sets that enables them to be smoothly joined to the 
best single-band PPTA data sets, giving data spans of 
up to 17.1 years. 

5 Future Prospects 

Figure [TTJ shows that the spectral index for the ex- 
pected stochastic GW background is generally steeper 
than the spectral index of other red processes that 
currently affect pulsar timing residuals. Furthermore, 
most of these red processes are different for different 
pulsars, whereas the variations due to the GW back- 
ground (and time-standard and ephemeris errors) are 
correlated. We can therefore be confident that, given 
sufficiently long data spans and a large enough sample 
of pulsars, PTA observations will reveal a GW back- 
ground unless current predictions of its strength at 
nanohertz frequencies are incorrect. 

Figure [12] shows the dependence of PTA sensitivity 
on the amplitude A s (Equation [TJ) of a stochastic GW 
background as a function of total data span T Q b s and 
number of pulsars in the PTA, -/V psr . These curves 
were obtained using Equation A12 of IVerbiest et al.l 
(20090 simulating an ideal PTA with N pBI pulsars 
distributed randomly on the sky, assuming weekly ToA 
sampling, 100 ns rms uncorrelated timing residuals, 
and no "red" noise contributions from sources other 
than the GW signal. 

The lowes t black line i s the P PTA "reference" case 
discussed by iJenet et al.l (|2005l ). Doubling the ob- 
serving time decreases the GW amplitude required for 
detection at a given significance level by a factor of 
approximately 4.5 ~ 2 13,/6 (cf., Equation [3] remem- 
bering that the significance is in terms of amplitude 

16 Notc that in lVerbiest et al.l <2009D Equation A8 the de- 
nominator should be squared and the number of spectral de- 
grees of freedom should be JVd.o.f. ~ l-^^obs/e where T ^ B 
is the data span and f c is the corner frequency where the 
white- noise power and gravitational- wave power are equal. 
The reasons for this change are: a) at least at higher signal 
levels, quadratic fitting does not remove a degree of free- 
dom from the fit, and b) the relevant bandwidth is that 
of the squared spectrum after Wiener filtering, that is, ap- 
proximately 0.7/c. 
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rather than power). The effect of doubling observing 
time on significance at a given amplitude greatly de- 
pends on the level of the GW amplitude; generally it 
is smaller for larger amplitudes because the GWs pass- 
ing over the pulsars contribute an uncorrelated noise 
term, thereby reducing the significance. Increasing 
the number of pulsars in the PTA improves the sig- 
nificance at a given GW amplitude by approximately 
[N psl (N psr — l)] 1 ^ 2 ~ N psT . Increasing the white noise 
rms amplitude by a factor of five, to 500 ns, has almost 
exactly the same effect as halving the data span, i.e., 
10 years with rms residual of 500 ns gives essentially 
the same sensitivity curve as 5 years with 100 ns rms 
residual. Given the sensitivity we can achieve with 
current instruments, the 10-year - 500-ns scenario is 
more realistic for the PPTA to achieve the necessary 
detection sensitivity. 

Improving rms timing residuals will also help, but 
this is difficult with present telescope facilities. In 
the Northern Hemisphere there are a number of radio 
telescopes with collecting area larger than that of the 
Parkes telescope which should give better results given 
similar instrumentation and similar observing cadence. 
For Parkes, a more sensitive broad-band receiver is cur- 
rently under development and this will give improved 
results. 

However, as Figure [12] makes clear, the greatest 
benefit is obtained by increasing the nu mber of pul- 
sars in the PTA. On-going searches (e.g.. Keith et akl 
l2012aT ) have already found MSPs which are suitable 
for PTAs and the best of these have already been in- 
cluded in existing PTAs. In time these will contribute 
fully to PTA sensitivity. This approximately linear in- 
crease in PTA sensitivity with N psT is the strongest 
motivation to combine the results from the three ex- 
isting (and possible future) PTA projects. The EPTA, 
NANOGrav and the PPTA are collaborating to form 
the In ternational Pulsar Timing Array (IPTA) (jHobbs et 
l2010al ). Although there is considerable overlap be- 
tween th e samples of the th ree main existing PTA 
projects, iHobbs et alJ ((2010a) list 37 pulsars that are 
regularly observed by at least one PTA. Of course, 
the data quality varies considerably across this sample, 
but there is no doubt that an IPTA data set will give 
superior results compared to the data set of any one 
participating PTA. E ven at the pessimist ic end of cur- 
rent predictions (e.g., ISesana et al.ll2008l ). IPTA data 
sets with spans of 10 years or more should give a sig- 
nificant detection of the GW background from binary 
black holes in distant galaxies. Failure to detect the 
GW background with these data sets will imply that 
current models for this background are fundamentally 
flawed. This would have important implications for 
current theories of black-hole formation and evolution 
in galaxies and for theories of galaxy growth through 
mergers. 

Even if a detection of GWs is achieved with IPTA 
data sets, it is clear that detailed study of the prop- 
erties of nanohertz GWs and their sources will have 
to w ait until the Squar e Kilometer Array (SKA) (see, 
e.g.. lCordes et al]|2004f ) is operational. The enormous 
increase in sensitivity provided by the SKA will allow 
detection and subsequent frequent timing observations 



of a large sample of MSPs. Ultimately a sample of 
200 PTA-quality MSPs should be possible, although it 
will take time to build up a significant data span. Once 
that is done, Figure[l2lshows that high-significance de- 
tetections will be possible even if the GW background 
is relatively weak, enabling accurate measurement of 
properties such as spe ctral index and possible orig in of 
the background (e.g., Ivan Haasteren et al.ll201lT ) and 
inves tigation of non- GR effects in the received signals 
(e.g. jLee et alj|2008h . 

Although Figure [T2] is computed for the case of 
detection of a stochastic GW background, similar con- 
siderations apply to other PTA objectives. Improved 
data sets should enable the detection and study of in- 
dividual GW sources, for e xample, bright black-hole 
binaries (Sesana et al. 2 009[) or burst si g nals from cos- 
mic strings (|Damour fc Vilenkinl 120051 ; ISanidas et alJ 
2012). Source locations with arc-minute precision or 
better should be measurable a nd other source prop- 
erties can be studied in detail (|Lee et al ]HqT3). Ex- 
tended data sets giv e improved pulsar-based timescales 
(jHobbs et al.l [2012? ) , enable improved determinations 
of the masses of solar-system planets and possibly even 
allow the detection of currently unknown solar-system 
objects such as TNOs. The future prospects of PTA 
research are indeed exciting. 
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Figure 1: Availability of the receivers and backend systems used for the PPTA project. The hatched area 
on the 10cm/50cm bar indicates when the 50cm was retuned to operate in the 700 - 770 MHz (40cm) band. 
The WBC bar is open when various instrumental problems affected the data quality. Significant intervals 
of overlap between operation of the various backend instruments allowed checks on instrument-dependent 
delays. 



Table 2: Backend systems used for the PPTA project 



System Bandwidth Dig. Bits Nr Channels Nr Bins P m ; 





(MHz) 








(ms) 


WBC 


1024 


2 


1024 


1024 


57 


PDFB1 


256 


8 


2048 


2048 


83 


PDFB2 


1024 


8 


1024 


2048 


4.1 


PDFB3 


2 x 1024 


8 


2048 


2048 


4.1 


PDFB4 


1024 


8 


2048 


2048 


4.1 


CPSR2 


2 x 64 


2 








APSR 


16 x 64 


2-8 
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Figure 2: Block diagram of the PDFB3/APSR system. The dashed line encloses components contained 
on the two main DFB boards. The 32-MHz synthesizer and samplers are on a separate board. In normal 
operation, the A and B inputs are used for the two polarisation channels from the receiver/IF system. 
The C and D channels may be used for the RFI adaptive filter reference input or independently for other 
signals. The L and H channels from the polyphase hlterbank refer to the lower and upper halves of the 
total bandwidth. Profiles from the pulsar binning memory are transferred to the control computer each 
DFB cycle. The APSR baseband outputs are output on two pairs of 10 Gb ethernet lines to switches which 
then distribute the signals amongst the 16 dual quad-core processors for quasi-real-time dedispersion and 
folding. The Control Computer has control links to most functional elements in the system, but most of 
these are omitted for clarity. 
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J1713+0747 1369 MHz J1713+0747 1369 MHz 




200 400 600 800 1000 




Pulse Phase Frequency (MHz) 

Figure 3: Total intensity (Stokes J) pulse profile displays formed using Psrchive routines for 
a 20cm PDFB4 observation of PSR J1713+0747. The upper-left plot is a false-colour image of 
the dedispersed pulse profile for each 1-min sub-integration, the upper-right plot shows the profile 
summed in time as a function of frequency across the band. The lower left plot is the mean pulse 
profile summed in time and frequency and the lower right plot shows the receiver bandpass for 
the two polarisations which are summed to form the total intensity. The upper plots show raw 
uncalibrated data whereas data for the mean profile plot has been band-pass and flux-density 
calibrated after excision of the few narrow RFI signals visible on the band-pass plot. A decrease 
in pulse intensity resulting from diffractive interstellar scintillation over the one-hour observation 
is visible in the upper-left plot whereas most of the frequency-dependent variations seen in the 
upper-right plot result from the instrumental bandpass. 
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Figure 4: Observed mean pulse profiles (red), fitted von Mises components (blue), the noise- free template 
obtained by summing the components (black) and, offset below the other profiles, the difference between 
the mean pulse profile and the template profile, for three of the PPTA pulsars at each of the three observing 
bands. The full pulse period is shown in each case. For PSR J0437— 4715, the mean pulse profiles are 
invariant interval; for the other two pulsars they are total intensity. 
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Figure 6: Distribution in celestial coordinates of pulsars suitable for pulsar timing array observations. All 
are radio-emitting MSPs (with P < 20 ms) in the Galactic disk except PSR J1824-2452A which is an 
MSP located in the globular cluster M28 (see text). The area of the plotted circle is inversely proportional 
to the pulsar period and the circles are filled for pulsars with mean flux density above 2 mJy. The dashed 
line is the northern declination limit of the Parkes radio telescope. Pulsars selected for the PPTA are 
marked with a star, red for the original 20 pulsars and mauve for the two pulsars recently added to the 
PPTA sample (see text). 



Table 3: The PPTA pulsars: basic parameters, observation times, flux densities and pulse widths 



PSR 
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DM 


Ph 


-^obs 


S700 
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S3 
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Wi 
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Mean 


Rms 


Mean 


Rms 


Mean 


Rms 
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(ms) 


(cm -3 pc) 


(d) 


(min) 


(mJy) ( 


mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(ms) 


(ms) 


J0437-4715 


5.757 


2.65 


5.74 


64 


406 


502 


149 


36 


32.3 


3.2 


0.1410 


0.0005 


J0613-0200 


3.062 


38.78 


1.20 


64 


7.2 


0.7 


2.3 


0.3 


0.42 


0.64 


0.462 


0.001 


J0711-6830 


5.491 


18.41 




64 


6.6 


4.5 


3.2 


5.7 


0.52 


0.26 


1.092 


0.005 


J1022+1001 


16.453 


10.25 
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64 


5.7 


2.5 


6.1 


5.1 


1.30 


0.30 


0.972 


0.005 


J1024-0719 


5.162 


6.49 




64 


5.4 


1.1 


1.5 


1.1 


0.37 


0.13 


0.521 


0.010 


J1045-4509 


7.474 


58.15 


4.08 


64 


9.2 


1.9 


2.7 


0.7 


0.48 


0.10 


0.840 


0.012 


J1600-3053 


3.598 


52.19 


14.34 


64 


3.4 


0.1 


2.5 


0.4 


0.77 


0.17 


0.094 


0.001 
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14.842 


38.05 


6.31 


64 


12.1 


2.9 


3.1 


0.9 


0.32 
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1.206 


0.003 
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62.41 


147.02 


64 
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0.314 


0.002 
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1.9 
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0.965 


0.004 


J1732-5049 


5.313 


56.84 


5.26 


64 


6.9 


2.7 


1.7 


0.3 


0.37 


0.06 


0.292 


0.003 


J1744-1134 


4.075 


3.14 




64 


7.8 


4.4 


3.1 


2.6 


0.71 


0.51 


0.137 


0.001 


J1824-2452A 


3.054 


119.86 




61 


10.6 


1.0 


2.0 


0.1 


0.33 


0.05 


0.972 


0.003 


J1857+0943 


5.362 


13.31 


12.33 


32 


10.5 


2.4 


5.0 


3.5 


1.01 


1.13 


0.518 


0.002 


J1909-3744 


2.947 


10.39 


1.53 


64 


6.1 


7.0 


2.1 


1.7 


0.77 


0.51 


0.0437 


0.0002 


J1939+2134 


1.558 


71.04 




32 


63 


19 


13.2 


5.0 


1.55 


0.72 


0.0382" 


0.0001° 


J2124-3358 


4.931 


4.62 




32 


12.3 


11.8 


3.6 


1.7 


0.44 


0.09 


0.524 


0.006 


J2129-5721 


3.726 


31.85 


6.63 


64 


5.0 


2.1 


1.1 


0.7 


0.11 


0.05 


0.262 


0.002 


J2145-0750 


16.052 


9.00 


6.84 


64 


16.4 


8.9 


8.9 


12.5 


1.38 


0.51 


0.337 


0.002 


a APSR data 



www.pu blish . csiro. au /journals/ pasa 




Figure 7: Timing residuals for the three PPTA bands, 50cm (red x), 20cm (black filled square) 
and 10cm (blue open circle) for four of the PPTA pulsars. Parameter files are from the 3-band 
solutions (see H4. 31 below) with the DM corrections and interband jumps set to zero and all other 
parameters held fixed. 
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Tabic 4: lOcm-band timing results for the PPTA pulsars 
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6.62 
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20 


0.37 
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0.66 
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56 


1.47 


ll c 
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1.16 


J1857+0943 
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a: Invariant interval, b: 2-year data span, c: 1.5-year data span, d: 6- month data span. 



Table 5: 20cm-band timing results for the PPTA pulsars 
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Data span 
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N iyr 


lyr rms res. 


xi 


Med. CTToA.lyr 




(MJD) 
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(ijs) 


J0437-4715" 


53431 - 55615 


998 


0.038 


95 


0.087 


7.48 


0.031 


J0613-0200 


53431 - 55620 


218 


0.88 


31 


0.81 


1.71 


0.63 


J0711-6830 


53431 - 55620 


212 


2.49 


36 


0.46 


0.99 


2.38 


J1022+1001 h 


53468 - 55618 


246 


1.17 


26 


0.66 


3.78 


0.84 


J1024-0719 b 


53431 - 55620 


175 


1.74 


28 


1.06 


1.48 


2.10 


J1045-4509 6 


53450 - 55620 


185 


2.13 


31 


1.75 


1.04 


1.81 


J1600-3053 fc 


53430 - 55598 


237 


0.50 


22 


0.44 


2.76 


0.29 


J1603-7202 h 


53430 - 55619 


168 


1.00 


27 


1.10 


2.42 


0.83 


J1643-1224 


53452 - 55598 


148 


0.69 


22 


0.88 


2.52 


0.60 


J1713+0747 


53452 - 55598 


198 


0.16 


23 


0.12 


5.86 


0.076 


J1730-2304 


53431 - 55598 


130 


1.20 


16 


0.98 


2.36 


0.87 


J1732-5049 


53724 - 55582 


102 


2.08 


19 c 


1.47 


0.78 


1.85 


J1744-1134 6 


53452 - 55599 


169 


0.38 


20 


0.17 


1.71 


0.30 


J1824-2452A 


53518 - 55620 


178 


0.48 


15 


0.57 


5.65 


0.32 


J1857+0943 6 


53431 - 55599 


121 


1.07 


18 


0.69 


1.33 


1.07 


J1909-3744 


53431 - 55620 


396 


0.13 


33 


0.097 


4.11 


0.086 


J1939+2134 


53450 - 55599 


139 


0.12 


gd 


0.096 


23.14 


0.035 


J2124-3358 


53431 - 55619 


184 


2.08 


25 


1.34 


1.63 


1.31 


J2129-5721 fc 


53476 - 55618 


182 


2.27 


28 


0.88 


0.91 


1.38 


J2145-0750 fc 


53431 - 55618 


215 


1.25 


29 


0.44 


1.43 


1.06 



a: Invariant interval, b: MEM calibration, c: 1.5-year data span, d: 6-month data span. 
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Table 6: 50cm-band timing results for the PPTA pulsars 



PSR 
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N lyr 
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xi 


Med. a T oA,iyr 




(MJD) 
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J0437-4715" 
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0.30 


72 


0.162 


2.83 


0.14 


J0613-0200 


53451 - 55619 


117 


0.53 


24 


0.73 


2.78 


0.45 


J0711-6830 


53431 - 55619 


127 


4.15 


26 


2.11 


1.38 


3.12 


J1022+1001 


53431 - 55619 


119 


1.65 


21 


1.17 


1.57 


1.40 


J1024-0719 


53431 - 55619 


75 


4.00 


20 


1.07 


0.58 


2.83 


J1045-4509 


53450 - 55619 


125 


3.00 


24 


3.37 


3.23 


2.11 
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53431 - 55619 


97 
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20 
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0.90 


1.43 
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17 
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0.99 


J1730-2304 
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84 
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16 


1.43 


1.70 


1.46 
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54 
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9 b 
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1.63 


4.16 
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53432 - 55620 


82 


0.66 


20 


0.65 


4.59 


0.37 


J1824-2452A 


53606 - 55583 


68 


1.16 


12 c 


2.64 


11.04 


1.16 


J1857+0943 


53482 - 55620 


74 


3.00 


17 


2.40 


1.48 


2.18 


J1909-3744 


53468 - 55619 
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0.30 


24 


0.39 


7.10 


0.24 


J1939+2134 


53484 - 55619 


79 
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10 d 


0.42 


42.15 


0.089 


J2124-3358 


53452 - 55619 


84 


3.64 


21 


1.27 
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2.33 


J2129-5721 


53686 - 55619 
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1.32 


22 


1.31 


2.60 


1.09 


J2145-0750 
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1.48 


20 


0.71 


0.63 


1.31 



a: Invariant interval, b: 2-year data span, c: 1.5-year data span, d: 6- month data span. 



Table 7: Timing Results for the PPTA pulsars 



PSR 


Band 


Corr. 


Data span 

(yr) 


N 


DMC Int. 

(d) 


Rms Res. 

M 


xi 


J0437-4715 


10cm 


IVI+DMC 


4.76 


18 


60 


0.075 


5.50 


J0613-0200 


20cm 


DMC 


6.00 


13 


700 


1.07 


1.76 


J0711-6830 


20cm 




6.00 


7 




0.89 


1.66 


J1022+1001 


20cm 


MEM 


5.89 


12 




1.72 


9.27 


J1024-0719 


20cm 


MEM 


6.00 


8 




1.13 


1.40 


J1045-4509 


20cm 


MEM+DMC 


5.94 


13 


70 


2.77 


1.80 


J1600-3053 


20cm 


MEM+DMC 


5.94 


13 


115 


0.68 


2.78 


J1603-7202 


20cm 


MEM 


6.00 


12 




2.14 


7.93 


J1643-1224 


20cm 




5.88 


14 




1.64 


5.46 


J1713+0747 


10cm 




5.71 


16 




0.31 


4.00 


J1730-2304 


20cm 


DMC 


5.94 


7 


300 


1.47 


2.90 


J1732-5049 


20cm 


DMC 


5.09 


12 


600 


2.22 


1.34 


J1744-1134 


20cm 


MEM+DMC 


5.88 


8 


1000 


0.32 


4.77 


J1824-2452A 


20cm 


DMC 


5.76 


7 


82.5 


2.44 


30.22 


J1857+0943 


20cm 


MEM+DMC 


5.94 


12 


425 


0.84 


1.16 


J1909-3744 


10cm 


DMC 


5.76 


17 


105 


0.133 


2.21 


J1939+2134 


20cm 


DMC 


5.88 


7 


50 


0.68 


141.63 


J2124-3358 


20cm 


DMC 


6.00 


8 


1800 


1.90 


1.38 


J2129-5721 


20cm 


MEM+DMC 


5.87 


12 


1500 


0.80 


1.00 


J2145-0750 


20cm 


MEM 


6.00 


14 




0.78 


3.18 
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Figure 8: Observed dispersion-measure (DM) variations for four PPTA pulsars 



Table 8: Pulse frequency sccond-timc-dcrivativcs for PPTA pulsars 



PSR 


v (3-band) 
(10- 28 s" 3 ) 


v (best) 
(10- 28 s- 3 ) 


PSR 


v (3-band) 
(10- 28 s- 3 ) 


v (best) 
(10- 28 s- 3 ) 


J0437-4715 


-1.9 ±0.9 


-5.3 ±0.6 


J1730- 


2304 


-11 ± 10 


-2 ±9 


J0613-0200 


1 ± 14 


8± 12 


J1732- 


5049 


55 ±57 


6 ±31 


J0711-6830 


-12 ±8 


-23 ±6 


J1744- 


1134 


5.7 ±3.8 


7.9 ±3.3 


J1022+1001 


10.6 ±2.9 


9.4 ±3.1 


J1824- 


2452A 


392 ± 48 


314 ±23 


J1024-0719 


-40 ±8 


-35 ±7 


J1857+ 


-0943 


-2 ±14 


4±8 


J1045-4509 


6 ±44 


20 ± 14 


J1909- 


3744 


7.8 ±3.5 


10.1 ±2.1 


J1600-3053 


30 ± 17 


39 ±6 


J19394 


2134 


231 ±22 


216 ±11 


J1603-7202 


-22 ±7 


-3 ±6 


J2124- 


3358 


-13 ±15 


5 ±14 


J1643-1224 


-100 ±37 


-53 ±15 


J2129- 


5721 


26 ± 14 


13 ±9 


J1713+0747 


12 ±4 


8±4 


J2145- 


0750 


11.9 ±3.3 


7.4 ± 1.7 
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Figure 9: Analytic (noise- free) timing profile templates for the "best" band for each of the 20 PPTA 
pulsars. The full pulse period is shown in each case and the vertical dashed line gives the template 
reference phase. 
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Figure 10: Final post-fit timing residuals for the PPTA pulsars for the band and corrections as listed in 
Table [7] The vertical extent of each subplot is adjusted to fit the data and its value is given below the 
pulsar name. The dashed line marks zero residual. 
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Figure 11: Power spectra of fluctuations in the timing residuals for the PPTA pulsars. The dashed line 
is the expected power spectrum in the timing residuals for a stochastic gravitational-wave background 
signal of amplitude A g — 1CP 15 . Note that the y-axis range is the same for all pulsars and covers nine 
orders of magnitude in power. Pulsars that have a i) with significance > 3er (Table [5]) are marked with * 
after the name. 
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GW Characteristic Amplitude 



Figure 12: Sensitivity of a PTA to a stochastic background of GWs as a function of total data span T Q b s , 
number of pulsars N in the PTA and assumming 100 ns rms timing residuals. Black lines are for N psi = 20 
and, T bs of 5 yr (unmarked), 10 yr (x) and 20 yr (o), respectively, red lines are similar for AT psr = 50 
and blue lines for N USI = 200. 
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Appendix A — Extension of the 
PPTA Data Sets 

Many of the PPTA pulsars had been observed at Parkes 
for some year s prior to the commencement of the PPT A 
project (eg., iToscano et al.|[l999l ; iHotan et~all 120061 1 . 
Long data spans are very beneficial to many pulsar 
timing studies and so it is useful to combine these ear- 
lier data with the PPTA observations. Here we de- 
scribe the "Extended PPTA" d ata sets that add ear- 
lier o bservations as described bv lVerbiest et afl (2008, 
120091 ) to the PPTA data sets. 

IVerbiest et al.l (|2008T l analysed Parkes observations 
of PSR J0437-4715 obtained over a 10-year period 
from 1996 to 2006 using two different receivers (H-OH 
and multibeam centre beam) operating in the 20cm 
band and with four different back-end systems oper- 
ating at different times. With these data, they ob- 
tained a precise distance for the pulsar, both with a 
direct measurement of the annual parallax and by mea- 
suring the apparent orbital period derivative resulting 
from the transverse motion of the system. The back- 
end systems us ed were the Fast Pu lsar Timing Ma- 
chine (FPTM) (ISandhu et al.l 11997ft in 19 96 - 1997, 
the S2 VLBI recorder ( van Stratenl [2001 ) in 1997 - 
1998, the Caltech -Parkes-Swinburne recorder CPSR 
(| van Stratenl l2003t ) for 1998 - 2002 and CPSR2 for 
2002 - 2006. 

Similarly, Verbi est et al.l (|2009l ) presented 20cm data 
for all 20 PPTA pulsars with data spans ranging from 
3.9 years to 14.2 years recorded between 1994 and 
2008. The FPTM was used f rom 1992 to 2001, a 5 12- 
channel analogue filterbank (|D'Amico et al.ll200ll ) in 
2002 and 2003 and CPSR2 from 2002 November to 
2008. For three pulsars, PSRs J1045-4509, J1909-3744 
and J1939+2134, data from the 50cm receiver was used 
to model the DM variations after this receiver became 
available in 2003 November. 

These data were analysed bv lVerbiest et al.l (|2008l . 

2009) using Tempo2 and fitting for the pulsar astro- 
metric, pulse frequency and binary parameters (if ap- 
propriate) simultaneously wi th arbitrary offsets be tween 
the different instruments. As lYardlev et all (|20 1 lT ) have 
shown, fitting for offsets (or jumps) subtracted most 
of the low-frequency power from the post-fit residu- 
als. We hav e therefore reanalysed the IVerbiest et all 
(2008, 2009) data sets to separately determine and fix 
as many of the offsets as possible. Where overlapping 
or near-overlapping (gap < 150 d) data from two differ- 
ent instruments were available, the offsets were mea- 
sured from short data spans (typically a few months) 
fitting for just v and the offset with all other param- 
eters held at the values from a fit to the entire data 
span for that pulsar. Offsets determined in this way 
were held fixed in subsequent analyses. For a few pul- 
sars, data gaps between some pairs of instruments were 
too large for this procedure to give reliable results and 
arbitrary jumps were retained. 

The lVerbiest et al.l ([2008, 200$ data sets were then 
combined with the three-band PPTA data sets de- 
scribed in the main text to form extended PPTA data 
sets with data spans of up to 17.1 years. These ex- 



tended data sets were then analys ed using Tempo2 
with the Cholesky decomposition (| Coles et al.l 1201 ll ) 
to properly handle the red noise component. The DM 
and DM offsets determined as described in i)4.3l for the 
PPTA data spans were incl u ded and hel d fixed. Since 
most of the IVerbiest et ail (|2008L 2009) observations 
were at frequencies in or close to the 20cm band, it 
was generally impossible to obtain DM variations from 
these data that were sufficiently precise to enable cor- 
rection for variable dispersion delays. The exception 
is PSR J1939+2134, where dual-band FPTM observa- 
tions (band centres about 1420 MHz and 1650 MHz) 
enabled a DM offset at MJD 50350 (1996 September) 
to be measured. This was held fixed in the subsequent 
Cholesky fit along with the other DM offsets. 

For all pulsars, parameters for which significant 
values could be obtained were included in the fit. Fi- 
nally, the remaining instrumental offsets were included. 
The use of the Cholesky method ensures that these off- 
sets (along with the other parameters) and their un- 
certainties are reliably determined despite the presence 
of the red noise. The first seven columns of Table [9] 
give the data spans in MJD and years, the total num- 
ber of ToAs (all bands), number of pulsar parameters 
fitted (N p ), the number of fixed instrumental offsets 
or jumps (Njo) and the number of instrumental offsets 
included in the final fit (Nji). The fixed instrumen- 
tal jumps include three interband jumps; the "best" 
band is designated as reference and the other three 
bands (including both 40cm and 50cm) are referred 
to it. Fitted jumps were needed for only seven pul- 
sars and for all of these except PSR J0437— 4715, the 
number required was three or less. PSR J0437— 4715 
required more offsets because more instruments were 
used for observations of this pulsar and the high preci- 
sion of the ToAs revealed a larger number of significant 
offsets. 

To form a "best" extended data set, the Verbies t et al.l 
(|2008l . 120091 ) data sets were combined with the best 
PPTA data set described in § 03] and Table [7] Pa- 
rameter files were created by copying the pulsar pa- 
rameters and jumps and DM corrections for the ex- 
tended three-band fit described above and holding all 
parameters except v and v fixed. Columns 8 - 10 of 
Table [5] give the number of ToAs in the best single- 
band extended data set, and the post-fit rms timing 
residual and reduced for each pulsar. Note that 
Cholesky decomposition was not used in this fit. Fig- 
ure [13] shows the post-fit residuals, clearly illustrating 
the better quality of the PPTA data sets compared 
to those obtained with the earlier systems. For several 
pulsars, e.g., PSRs J0437-4715 and J1024-0719, they 
also reveal red noise that was obvious in neither the 
6-year PPTA data set s (Figure | 10|l. nor in the po st-fit 
residuals presented bv lVerbiest et al.1 i|2008l , 120091 ). 

To quantify these in the same way as for the PPTA 
data sets (Table [8} , v was additionally fitted, giving 
the values in the last column of Table O In contrast 
to the PPTA data sets where ten of the 20 pulsars had 
v values significant at the 3<r level, all but three of the 
pulsars have significant v values for the extended data 
set. This is a simple illustration of the importance of 
long data spans for characterising red-noise processes. 
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PSR J1024— 0719 is a particularly striking example: 
for the PPTA data set this had a v with a significance 
of 5<J, whereas for the extended data set the value is 
consistent but the significance has increased to 48<r. 
It is interesting to note that, based on the PPTA data 
set ( i)4,5[) . this pulsar appears to have a relatively steep 
red timing noise spectrum. 

PSRs J1824-2452A and J1939+2134 stand out with 
positive and relatively large v values, but it is no- 
table that 13 of the 17 significant v in Table [S] are 
also positive. A small glitch was observed for P SRs 
J1824-2452A in 2001 ijcoenard fc Backer! l2004h . It 
is well known that large positive va lues of v are ob- 
served after glitch es in young pulsars l|Lvne et al ]|2000l : 
lYuan et al.ll2010T ). These are attributed to the dynam- 
ics of superfluid vortices in the interior of the neu- 
tron star as they re-establish equilibrium with th e rigid 
crust following a glitch fe.g.. lAlpar et af]|l993T ). This 
suggests the intriguing possibility that the observed 
positive v values are attributable to post-glitch recov- 
eries, in all cases, except possibly for PSR J1824— 2452A, 
from unseen earlier glitches. It is worth noting that the 
v values seen for the PPTA MSPs are about four or- 
ders of magnitude smaller than those seen in young 
pulsars. 

Comparison of the observed rms timing residuals 
fo r the extended data s ets wit h the "CP+MSP" model 
of IShannon fc Cordesl (2010) gives similar results to 
those for the PPTA data sets discussed in H4.5I Com- 
pared to the PPTA data sets, the (, values for most 
pulsars are a factor of 2 - 8 smaller, indicating the di- 
minished relative contribution of the white noise with 
the longer data spans. For PSR J0437-4715, C ~ 0.15 
showing that the model substantially over-predicts the 
amount of timing noise with this longer data span. For 
PSRs J1909-3744 and J1939+2134, £ values are 1.18 
and 0.49 respectively, similar to those for the shorter 
data spans; for PSR J1824— 2452A there are no new 
data in the extended data set. 

For several pulsars, most notably PSRs J0437— 4715 
and J1713+07 1 3, there is significant red noise in the 
IVerbiest et al.l (2008, 200^) data sets al one which was 
not present in the results presented bv I Verb icst et al. 
(2008, 120091 ). This is mostly attributable to the fix- 
ing of the instrumental offsets. However, it is likely 
that the inability to accurately model DM variations 
for most of the early data also contributes significantly 
to the observed residuals. For example, the deviations 
observed in the early data for PSR J 1045— 4509 are al- 
most certainly due to uncorrected DM variations (cf., 
iKeith et alJl2012bh . It is possible that observations 
from other observatories will help to determine pre- 
PPTA DM variations. However it is generally difficult 
if not impossible to establish absolute ToAs for archival 
data with sufficient precision and so these data are 
likely to be of limited use in this respect. Only fu- 
ture carefully calibrated multi-band observations will 
establish for sure if DM variations are a significant 
contributor to the observed red noise on decade-long 
timescales. 



www.pu blish . csiro. au /journals/ pasa 



35 



Tabic 9: The Extended PPTA data sets 



PSR 




Data Span 


NtoA,3B 


N P 


Njo 


Nji 


NtoA 


Rms Res. 


xi 






(MJD) 


(yr) 












M 




(10- 28 s" 3 ) 


J0437-4715 


50190 


- 55619 


14.86 


5055 


15 


7 


9 


3508 


0.21 


7.17 


1.26 ±0.03 


J0613-0200 


51526 


- 55619 


11.21 


629 


12 


6 





341 


1.11 


1.24 


7.2 ±2.1 


J0711-6830 


49373 


- 55620 


17.10 


555 


6 


10 


1 


319 


1.54 


1.54 


-0.8 ±0.7 


J1022+1001 


52649 


- 55619 


8.13 


624 


12 


7 





378 


1.82 


8.14 


-0.4 ±1.6 


J1024-0719 


50117 


- 55620 


15.07 


493 


6 


10 





309 


4.38 


12.57 


-38.6 ±0.8 


J1045-4509 


49405 


- 55620 


17.02 


635 


10 


10 





393 


5.05 


3.18 


9.3 ± 1.2 


J1600-3053 


52301 


- 55598 


9.03 


704 


12 


7 


1 


503 


0.98 


1.21 


8.6 ±2.2 


J1603-7202 


50026 


- 55619 


15.31 


483 


12 


7 


3 


290 


2.12 


3.08 


1.2 ±0.4 


J1643-1224 


49421 


- 55598 


16.91 


477 


11 


7 


3 


288 


2.30 


5.90 


6.0 ± 1.0 


J1713+0747 


49421 


- 55619 


16.97 


612 


15 


10 





334 


0.46 


7.75 


-2.60 ±0.16 


J1730-2304 


49421 


- 55598 


16.91 


390 


7 


10 





223 


2.59 


3.25 


-0.8 ±0.8 


J1732-5049 


52647 


- 55582 


8.04 


244 


11 


9 





149 


2.47 


1.17 


28 ± 7 


J1744-1134 


49729 


- 55599 


16.07 


534 


7 


9 


3 


368 


0.65 


3.27 


1.9 ±0.3 


J1824-2452A 


53518 


- 55620 


5.75 


302 


6 


3 





178 


2.02 


14.50 


241 ± 22 


J1857+0943 


53086 


- 55599 


6.88 


291 


15 


7 





152 


0.96 


1.18 


7.1 ±5.2 


J1909-3744 


52618 


- 55619 


8.22 


1245 


14 


7 





724 


0.19 


5.06 


3.54 ±0.44 


J1939+2134 


49956 


- 55599 


15.45 


386 


7 


9 


3 


237 


4.27 


3664 


127.8 ± 1.4 


J2124-3358 


49489 


- 55619 


16.78 


652 


7 


11 





473 


2.92 


1.85 


-6.1 ±0.8 


J2129-5721 


49987 


- 55619 


15.42 


448 


11 


10 





285 


1.41 


2.21 


6.3 ±0.8 


J2145-0750 


49517 


- 55618 


16.70 


972 


13 


10 





696 


1.06 


2.81 


-1.38 ±0.12 



36 



Publications of the Astronomical Society of Australia 



J0437-4715 

3.65yUS 




J0613-0200 
11.77^s 



J071 1-6830 
40.40Ats 



J1022+1001 
23.71jU.s 



J1024-0719 
55.20jUS 




J1045-4509 
114.83/as 



J1600-3053 
12.34jits 



J1603-7202 
28.34jLts 



Li 



J1643-1224 
45.09jus 



J1713+0747 
7.98^s 



ij # - 



J1730-2304 
23.22/xs 



sit 



J1732-5049 
28.15/^s 



J1744-1134 
9.95/xs 




J1824-2452A 
17.59/xs 



J1857+0943 
13.71jLts 



J1909-3744 
7.16yus 



J1939+2134 
42.02/xs 




J2124-3358 
75.09,us 



J2129-5721 
36.11jus 



I I I I I I I I I I 



J2145-0750 
45.97jLts 



1995 



2000 2005 
Year 



2010 



Figure 13: Post-fit timing residuals for the extended PPTA data sets. The vertical extent of each subplot 
is adjusted to fit the data and its value is given below the pulsar name. The dashed line marks zero 
residual. 



